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IN  case  the  main  conduit  leading  to  a  turbine  is 
comparatively  long,  beginning  at  the  intake  on  a 
comparatively  light  slope  (canal,  tunnel  or  low- 
pressure  pipe)  and  concluding  in  a  short  steep  slope 
leading  directly  to  the  turbines,  (high-pressure  pipe  or 
penstock),  it  is  common  practice  to  construct  at  the  junc¬ 
tion  an  open  tank,  called  the  surge  tank,  by  which  the 
main  conduit  is  separated  into  two  parts,  in  which  with 
a  constant  discharge  of  the  turbines  the  flow  in  the  con¬ 
duit  is  constant,  and  with  a  variable  turbine  discharge  a 
variable  conduit  flow  results.  In  the  first  case,  the  water 
surface  in  the  surge  tank  is  lower  than  the  water  surface 
at  the  intake  by  an  amount  dependent  upon  the  flow  in 
the  low-pressure  pipe.  This  difference  of  head  is  deter¬ 
mined  by  the  friction  head  in  the  main  conduit.  In  an 
increment  of  time  just  as  much  water  flows  into  the  surge 


I.  Introduction. 

The  investigation  is  developed  (see  Fig.  i)  under  the 
following  additional  assumptions  : 

(1)  The  intake  is  provided  with  a  spillway  whose  di¬ 
mensions  are  such  that  the  elevation  n  —  n  in  the  forebay 
may  be  considered  as  constant  during  the  period  under 
investigation. 

(2)  The  sectional  area  of  the  main  conduit  is  constant. 

(3)  The  volume  of  the  conduit,  compared  to  the 
volume  of  the  surge  tank  and  compared  to  that  part  of 
the  volume  of  the  forebay  which  is  affected,  is  so  large 
that  the  influence  of  both  of  these  masses  of  water  toward 
decreasing  the  flow  may  be  neglected. 

(4)  Elastic  and  temperature  conditions  are  neglected. 

In  the  derivations  of  the  formulae,  the  following  ab¬ 
breviations  are  used  : 


tank  as  flows  out  of  it  through  the  penstocks  to  the 
turbines. 

In  the  second  case,  the  inertia  of  the  moving  mass 
in  the  main  conduit  prevents  this  equality  of  inflow  and 
outflow.  The  water  surface  in  the  surge  tank  has  variable 
tO  heights,  that  is,  it  rises  or  falls  above  and  below  the  ele- 
vation  due  to  the  steady  flow.  The  extent  of  this  fluctua¬ 
tion  of  the  water  surface  depends  upon  the  dimensions  of 
the  main  conduit,  the  amount  of  flow,  and  also  on  the 
size  and  form  of  the  surge  tank.  In  this  case,  if  an  over¬ 
flow  or  an  excessive  lowering  or  too  large  fluctuations  of 
the  water  surface  are  to  be  prevented,  the  surge  tank 
must  be  dimensioned  according  to  the  area  and  length 
of  the  main  conduit  and  according  to  the  inflow  and  out¬ 
flow.  In  this  article  the  problems  leading  to  the  deter¬ 
mination  of  the  dimensions  will  be  discussed,  and  the 
methods,  partly  analytical  and  partly  graphical,  involved 
in  the  solution  of  these  problems,  will  be  developed  upon 
the  assumption  that  the  main  conduit  from  the  intake  to 
the  surge  tank  consists  of  a  pipe  or  tunnel  running  full 
under  pressure. 

*  Reprinted  from  The  Canadian  Engineer,  Vol.  27. 


L  =  Length  of  main  conduit  in  feet. 

a  =  Sectional  area  of  main  conduit  in  square  feet. 

p  =  Wetted  perimeter  of  main  conduit  in  feet. 

v  =  Velocity  of  water  in  main  conduit  in  feet  per  second 
at  the  time  “t”. 

v,  =  Normal  velocity  of  water  in  the  main  conduit  in  feet 
per  second,  during  the  period  of  steady  flow. 

Vi  =  Initial  velocity  in  the  main  conduit  in  feet  per 
second,  at  the  time  t  =  o. 

(v,  Vi  and  v2  are  average  values  and  are  assumed  as 
constant  throughout  the  entire  length  of  the  main 
conduit). 

H ,  =  Vertical  distance  from  the  water  surface  n  —  n  in 
the  forebay  to  the  centre  of  gravity  of  the  entrance 
to  the  main  conduit. 

Hi  =  Vertical  distance  from  the  water  surface  n  —  n  in 
the  forebay  1o  the  centre  of  gravity  of  the  entrance 
to  the  surge  tank. 

Q ,  o,,  Qi  =  Volumes  of  water,  corresponding  to  time  t, 
the  time  during  steady  flow,  and  for  lime  t  —  a 
(beginning). 
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h,  h,  and  h2  =  Friction  heads  corresponding  to  the  pipe 
dimensions  and  velocities  v,  v,  and  v2. 
z  ~  Vertical  distance  from  the  elevation  n  — n,  of  the 
water  surface  in  the  surge  tank  at  the  time  t,  taken 
positive  above  n  —  n  and  negative  below  n  — n. 

A  =  Sectional  area  of  surge  tank  in  square  feet  for  the 
elevation  determined  by  z,  so  that  in  general  A  is  a 
function  of  z. 

s  —  Velocity  of  the  water  surface  elevation  in  the  surge 
tank,  positive  when  rising,  negative  when  falling, 
as  an  average  value  assumed  constant  in  the 
section  A. 

q  =  Volume  of  discharge  through  the  penstock  at  the 
time  t,  in  cubic  feet  per  second. 

q 

c  =  —  =  the  discharge  velocity  of  the  volume  q,  in  feet 
A 

per  second,  with  respect  to  A. 

More  notations  will  be  introduced  during  the  course 
of  the  investigations. 

II.  Derivation  of  Principal  Equations. — Referring  to 
Fig.  2,  the  distance  between  two  adjacent  cross-sections 

w  .a. dl 

of  the  main  conduit  =  dl,  and  therefore  -  =  mass 

g 

of  the  water  between  the  two  sections  ( w  =  weight  of 
cubic  unit  of  water,  g  =  acceleration  due  to  gravity). 

At  the  left  side  cross-section  at  the  time  t,  there 
exists  the  pressure  p  in  pounds  per  square  foot,  an 
average  value,  assumed  constant  for  the  entire  cross- 
section.  In  the  right,  the  pressure  is  pr  at  the  same 
time  t.  The  whole  of  p'  is  in  general  different  from 
that  of  p,  by  the  amount  dp.  The  pressure  P  is  a  func¬ 
tion  of  the  location  of  the  cross-section,  that  is  to  say, 
depends  upon  l ;  (l  =  the  distance  of  the  left  side  cross- 
section  from  the  entrance  of  the  main  conduit)  and  de¬ 
pends  further  upon  the  time  t  (the  flow  varies  with 
the  time). 

Sp  8p 

dp  =  —  dl  +  — .dt  -  -  (i) 

SI  H 

Since  p  and  p'  occur  at  the  same  time,  p1  is  different 
from  p  by  the  difference  due  to  the  distance  between  the 
cross-sections,  and  therefore  in  the  foregoing  formula  the 

8  p 

differential  dt  of  the  time  =  o.  Therefore,  dp  =  — .dl  (2). 

H 

The  following  forces  act  upon  an  element  of  unit  mass,  in 
the  direction  of  flow,  that  is,  in  the  direction  of  v : 

1st.  The  weight  component  Pi  =  w.a.dl  sin «  (3) 
where  “  =  the  inclination  of  the  axis  of  the  main  conduit 
from  the  horizontal,  dl.  sin  a  —  cl  H  =  vertical  distance 
between  the  centres  of  gravity  of  the  sections  a  b  and  a,  b,. 

Pi  =  w.a.d  H  -  -  (4) 

2nd.  The  difference  between  the  reactions  due  to  the 

8  p 

pressure  p  and  p  +  — .dl  which  is 
SI 

8  p  8  p 

P2  —  p.a  —  (p  +  — .dl)  a  =  ■ —  a  — dl  -  (15) 

8 1  8  Z 

Contrary  to  the  direction  of  motion  the  friction  exerts 
a  force.  Let  K  be  the  amount  of  the  friction  in  terms  of 


the  velocity  v  and  we  get  P  =  —  wa .dl.K  (6)  where  K  re¬ 
presents  a  pure  number  with  respect  to  its  dimension. 

According  to  the  general  fundamental  law :  mass  x 

dv 

acceleration  =  acting  force,  it  results  (as  —  is  the  ac- 

dt 

celeration  with  respect  to  the  velocity  v  at  the  time  t)  that, 
dv 

m.  —  =  P.  +  P2  —  P  -  -  -  (7) 

dt 

w.a.dl  8  p 

-  =  w.a.dH  —  a  —  dl  —  w. a  dl.K  (8) 

g  *1 

dl  dv  1  5  p 

—  .  —  =  dH - dl  —  K  dl  -  -  (9) 

g  dt  w  H 

The  velocity  v  has  the  same  value  at  the  time  t  in  the 
whole  length  of  the  main  conduit.  The  case  is  therefore 
dv 

the  same  for  —  and  K.  Therefore,  if  we  integrate  be- 
dt 

tween  the  limiting  values  l  =  o  and  l  —  L  or  H  =  H, 
and  H  =  H2  relative  to  a  motion  from  the  intake  to  the 
main  conduit,  we  get  the  following:: 

L  dv  1  P  l  t>p 

~.—  =H2  —  H  1 - /  —  dl  —  KL  (10) 

g  dt  w  J  0  8/ 

A  hP 

In  the  integral  /  —  .  dl,  p  is,  as  demonstrated, 

J  0  82 

a  function  of  t  and  l.  But  as  the  integration  relates  to 
the  condition  at  a  certain  time  f,  t  is  to  be  considered  as 
a  constant  and  therefore 

L  Sp 

—  dl  —  p2  —  p>  -  -  (11) 

SI 

where  p>  =  pressure  at  the  entrance  to  the  surge  tank 
p,  at  the  entrance  to  the  main  conduit,  KL  is  nothing 
other  than  the  friction  head  h  for  the  entire  main  conduit 
at  the  time  t.  Now,  we  may  say  that 

P  2  pa  pi  Pa 

—  =  H2  +  z  +  — ;  —  =  Hi  +  —  (12) 

w  w  w  w 

p° 

in  which  —  equals  the  water  column  equivalent  to  the 
w 

atmospheric  pressure,  and  the  following  equation  results : 
L  dv 

—  .  —  +  z  +  b=  o  -  -  ( 1 3) 

g  dt 

Having  assumed  a  flow  in  the  main  conduit  from  the 
intake  to  the  surge  tank,  we  consider  now  the  flow  as  re¬ 
versed  (from  the  surge  tank  to  the  intake)  and  keep  the 
direction  for  the  measurement  of  the  length  l  the  same, 
then  we  must  consider  that  the  friction  now  acts  in  the 
direction  of  an  increase  of  l ,  contrary  to  the  first  instance. 
In  this  case  we  get  the  equation 

L  dv 

— - 1-  z  —  h  =  o  -  -  (14) 

g  dt 


CASE  A:  SUDDEN  SHUT-DOWN. 
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For  any  case  we  may  combine  both  equations  into  a 
fundamental  formula: 


If  we  introduce  as  an  abbreviation: 


L  dv 

—  . - h  z  ±  h  =  o 

g  dt 


H5) 


in  which  the  plus  refers  to  the  velocity  from  intake  to 
surge  tank  and  the  minus  sign  for  the  reverse  movement. 
If  we  now  introduce  for  h  a  function,  the  value  of  which 
has  the  same  sign  as  v,  then  the  double  sign  may  be 
dropped  and  the  whole  movement  may  be  represented  by 
the  equation : 

L  dv 

—  . - p  z  +  h  =  o  -  -  (16) 

g  dt 

A  second  fundamental  equation  may  be  developed 
from  the  condition  of  continuous  flow,  that  is  to  say,  the 
volume  of  water  which  flows  to  the  surge  tank  in  the  time 
element  dt  must  be  equal  to  the  sum  of  the  changes  of 
volume  in  the  surge  tank  and  in  the  volume  flowing  out 
of  the  tank  at  the  same  time 


v.a.dt  =  s  A  dt  +  q.dt. 

a.  v  =  As  +  q  =  T(s  +  c) 


(i7) 


The  questions  which  are  of  principal  interest  are: 

1  st.  How  is  the  movement  of  the  water  surface  in 
the  surge  tank  with  regard  to  time  affected  by  given 
dimensions  of  the  main  conduit  and  the  surge  tank  and 
given  conditions  as  to  overflow? 

2nd.  What  dimensions  must  the  surge  tank  have  to 
agree  with  given  conditions  as  to  the  main  conduit  and 
outflow,  in  order  that  the  fluctuation  of  water  level  does 
not  exceed  certain  amounts,  determined  by  local  condi¬ 
tions?  The  following  cases  will  be  investigated: 

(a)  Sudden  partial  or  entire  cessation  of  the  outflow. 

(b)  Sudden  starting  of  the  outflow. 

(c)  Gradual  cessation,  or  starting,  variable  outflow¬ 
ing  conditions;  and 

(d)  Influence  of  a  spillway  built  in  the  surge  tank. 

The  formulae  will  be  developed  in  the  first  case  for  a 

constant  sectional  area  of  surge  tank,  and  under  the 
assumption  that  h  is  proportional  to  v,  i.e.,  h  —  n.v. 

It  will  be  shown  that  it  is  possible  to  investigate  all 
cases  according  to  a  uniform  method,  analytically  and 
graphically,  with  the  assistance  of  the  well-known  theory 
of  damped  and  forced  undulations.  The  inaccuracies 
which  result  because  h  is  proportional  to  v 2  and  the  cor¬ 
rections  which  must  be  applied  to  the  results  of  the  first 
method  will  be  shown  in  the  final  study.  With  the  simpli¬ 
fying  assumptions  already  mentioned,  we  get  from 
equation  (17) 


(18) 


This  applied  in  equation  (16)  and  the  whole  equation 
divided  by 


dv 

A  ds 

dc 

—  = 

-  —  (— 

+  -1 

dt 

a  dt 

dt 

L 

A 

ds 

dc 

z 

T2  = 

—  . 

— 

(iq),  we 

get 

—  + 

■ — 

+  — 

g 

a 

dt 

dt 

T2 

ds 

n . 

A 

z 

n  A 

dc 

— 

+  — 

.  s  - 1 - 

+ 

.  c 

+  — 

dt 

T2 

.  a 

T2 

T2a 

dt 

s  = 


To 

T\  a 

n  .  A 

that 

dz 

ds 

d2z 

dt 

dt 

dt2 

(21) 


then  follows  the  principal  equation  : 

d2z  1  dz  2  c  dc 

-  +  —  .  —  +  —  +  —  4-  —  =  o 

dt 2  To  dt  T2  To  dt 

~L 


22 


(23) 


The  values  n  and  T  = 


*/- 


T2  a  L 

and  To  =  —  .  —  =  - 

n  A  n.g 


(24) 


are  times  with  regard  to  their  dimensions.  Starting  with 
this  principal  equation,  we  may  investigate  the  different 
cases  as  follows : 

III.  Special  Cases. 

Case  A.- — Sudden  Shut-down. 

Preceding  a  shut-down,  Q ,  cubic  feet  per  second 
flows  out  of  the  surge  tank.  During  the  normal  condition 

0. 

in  the  main  conduit,  the  velocity  —  Vi  =  —  wherefore 

a 

q  —  Q 1.  The  water  surface  in  the  surge  tank  is  h ,  feet 
lower  than  the  static  level  n  —  n.  The  time  t  is  measured 
from  the  beginning  of  the  shut-down,  therefore  from 
t  =  o,  q  becomes  eq  if  6  is  the  proportion  of  the  steady 
flow  subsequent  to  the  shut-down  in  relation  to  the  flow 
preceding  the  shut-down. 

After  the  sudden  shut-down,  the  following  phenomena 
occur  in  the  surge  tank  :  The  water  surface  rises  with 
variable  velocity  until  it  reaches  a  maximum  height. 
When  the  highest  elevation  is  reached,  the  reverse  move¬ 
ment  occurs.  The  velocity  increases  as  the  water  level 
recedes,  then  decreases  until  the  lowest  level  is  reached, 
after  which  an  ascending  movement  occurs  but  to  a  some¬ 
what  less  height  than  before,  and  so  on,  until  the  normal 
conditions  with  the  constant  flow  *.Q ,  become  established. 
The  movement  of  the  water  level  belongs  in  the  category 
of  the  damped  oscillations. 

(1)  Analytical  Investigation. — The  formula  (23)  for 
this  case,  under  the  condition  that 

q  6  Qi  dc 

c  =  —  =  -  =  e  c,  =  a  constant  ( —  =  o)  (25) 

A  A  dt 

becomes 

d2z  1  dz  z  Ci 

- 1 - 1 - +  e —  =  0  -  (26) 

dt 2  To  dt  T2  To 


Since  hi 


=  n.Vi  =  nc  1 


«  Ci 


T2 

T 


=  Ci  — ,  it  follows  that 


To 


« h, 


T2 


(27) 
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dz 

If  we  introduce  z  —  y  —  *hi,  therefore  — 

dt 

d2z  d2y 


dt 2  dt 2 


dy 

dt 

(28) 


from  (26)  follows  a  differential  equation  of  the  second 
order 

d2y  1  dy  y 

-  +  —  .  —  +  —  =  0  -  (29) 

dt 2  To  dt  T2 

This  is  a  linear  differential  equation  of  the  second 
order  with  constant  coefficients.  For  the  solution  we  say : 


x  1 

r2  +  —  r  +  —  =  0 

To  T2 

The  solution  of  this  quadratic  equation  will  give  us 
the  values  of  the  differential  equation  : 


We  know  that  according  to  the  theory  of  the  linear 
differential  equations  three  solutions  may  be  found. 

If  the  radical  is  positive  or  zero,  the  solution  of  the 
differential  equation  represents  a  non-periodical  function, 
that  is : 

t  t 


2  To  2  To 

y  =  Ri .  &  -}  R2 1  c 


1  I  I 

if  —  = - =  o  (30 a) 

T2  T2  (2  To)2 


t  t  t 


y 


+  — 

T, 

(I?i  e  +  R2  b 


T,  2  To 

)« 


But, 


1  1  1 

T,3  T3  (2  To)3 


negative  (30b) 


y  =  R  e 


2  To 

.  sin  (0  +  t/T  1) 


1 


if  — 
T2 


1  1 

T3  (2  To)2 


positive  (30c) 


represents  a  damped  harmonic. 

These  three  conditions,  due  to  the  different  values  of 
the  radical  of  the  quadratic  equation  can  be  written  also, 
substituting  for  T  and  T0  the  values  of  equations  24: 

1  A  4  7. 

—  becomes  zero  if  T  =  2  To  or  —  =  -  (31a) 

T2  a  n2g 

1  A  4  L 

—  becomes  negative  if  T  >  2T0  or  —  >  -  (31b) 

T2  a  n'g 


1  A  4  L 

—  becomes  positive  if  T  <  2T0  or  —  <  -  (31c) 

T,3  a  n2  g 

As  mentioned  before,  equations  30a  and  30b,  with 
the  conditions  of  31a  and  31b,  represent  non-periodical 
function,  that  is : 

Form  1  is  the  expression  for  a  damped  oscillation ; 
Forms  2  and  3  represent  non-periodical  movements, 
i.e.,  a  transition  from  one  quiescent  level  to  another  with¬ 
out  any  oscillations. 

As  will  be  seen  later  from  an  example,  n  has  in  most 
cases  a  value  which  lies  between  2  and  1  seconds ;  there¬ 
fore,  the  condition  for  a  non-periodic  water  level  fluctua¬ 
tion  is : 

A  L 

—  is  equal  to  or  greater  than  — . 
a  35 

Now,  for  L  —  x  miles 

A  should  be  equal  to  or  greater  than  150.X.U. 

This  case  may  well  occur  when  a  pond  is  used  as  the 
surge  tank.  With  artificially  constructed  surge  tanks 
A 

—  is  considerably  smaller;  the  following  investigations 
a 

are  therefore  limited  to  the  first  form  of  oscillations. 
As  z  =  y  —  eb,,  equation  (30)  may  be  written 

t 

2  To 

z  =  —  ebi  +  Re  sin  (3  +  t/Tx)  -  (32) 

dz 

and  the  differentiation  with  respect  to  t  gives  (since —  =  s) 


2  To  1  2  To  ] 

s  —  Re  -  - cos(0  +  f/T,)  —  sin  (£  + i/T,)  M33) 


2  To  T, 


2  To 


If  we  substitute  tgy  =  -  and  consider  that 

T, 


we  get : 


1  1  1 


T,  T3  4  To3 
t 


(34) 


R  2  To 

s  =  —  .  e  sin  (7  —  0  —  t/T,)  (35) 

T 

and  the  integration  constants  R  and  0  are  determined 
from  the  initial  conditions,  i.e.,  from  the  location  and  the 
condition  of  movement  of  the  water  level  in  the  surge 
tank  at  the  time  t  =  o. 

It  will  be  noted  from  the  description  of  the  phenomena 
for  that  case  that  for  the  shut-down  and  in  the  very  mo¬ 
ment  of  it,  the  water  surface  is  at  the  distance  h,  under 
the  elevation  n  —  n.  Therefore,  for  t  =  o,  z  —  z0  =  —  h,. 
The  initial  value  for  5  =  So  at  the  time  t  =  o,  must  be 
assumed  as 

Q>  —  €0. 

So  =  -  =  (1  —  €)  C,  -  (36) 

A 

This  assumption  does  not  only  consider  a  sudden 
shut-down,  but  assumes  also  a  sudden  beginning  of  a 
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uniformly  distributed  velocity  in  the  water  surface  of 
surge  tank.  This  assumption  is  true  for  the  limiting  cases 
because  a  short  duration  of  the  shut-down  does  not  ap¬ 
preciably  affect  the  results,  as  we  shall  see  from  the  in¬ 
vestigation  of  the  influence  of  a  shut-down  of  short 
duration.  Here  follow  the  equations  for  the  determina¬ 
tion  of  the  integration  constant  when  /  =  o  : 

R  sin  0  =  _  (i  —  «)  hx  -  -  (37) 

R  sin  (7  —  0)=  +  (1  —  •)  CiT  -  (38) 

The  latter  equation  may  be  transformed  to: 

T«  1 

R  cos  0  =  (1  —  «)  [ - ]  hi  T,  (39) 

T  2  To 


By  introducing  the  values  in  the  bracketed  term,  we 
find  that  this  term  is  proportional  to  the  difference 
a  n2  g 

( - )  -  -  (40) 

A  2  L 


A 

The  difference  is  therefore  positive  if  —  is  less  than 

a 


2  L 


n  g 


,  i.e. ,  with  the  same  assumptions  for  n  and  g  as  before : 
A  <  17.5 .x.a 


This  has  influence  on  the  determination  of  the  value 
of  0,  as  in  the  case  mentioned  sin  0  is  negative  and  with 
A  2  L 

—  <  -  the  cos  0  becomes  positive,  i.e.,  0  must  lie  in 

a  n*g 

4  L  A  2  L 

the  fourth  quadrant.  If -  >  —  >  - ,  the  bracketed 

n  g  a  n  g 

expression  and  .’.  cos  0  must  be  negative,  0  lies  in  the 
third  quadrant.  For  the  last  formula,  we  may  find  for 
the  determination  of  R  and  0  the  equations : 

Ti  To 

R  =  (1 — 6)  ht -  -  -  (41) 

T 


1 

tg0  - - -  (42) 

T./T,  —  %  T,/T„ 


This  value  used  for  equation  32  gives  the  following 
maximum  values : 


z  max,  =  —  +  Re 


T, 

- (7  —  0) 

2  To 


sin  7 


Ti 

- -  (7  - 0  +  2*) 

2  To 

z  max,  =  —  6/i,  +  Re  sin  7,  etc.  (44) 

and  the  following  minimum  values: 

T, 

- (7  —  0  +  «■) 

2  To 

z  min,  =  —  «ft,  —  Re  sin  7 


T, 

—  - -  (7  —  0  +  3  ») 

2  T. 

z  min,  =  —  —  Re  sin  7,  etc.  (45) 

Due  to  z  =  y  —  «Ji, 
y  max,  y  max, 


y  max,  y  max. 


T, 

- .  2* 

2  To  y  min,  .  y  min, 

=  e  =  - - —  =  - —  = - t46) 

y  min,  y  min. 

The  amplitude  of  this  oscillatory  motion  is  decreasing. 

Iff  =  infin.  z  becomes  —  *h, 

5  becomes  zero 

(2)  Graphical  Demonstration. — For  the  graphical 
demonstration  of  the  velocity,  a  very  simple  construction 
may  be  used.  We  may  divide  the  equation  32  into  two 
equations  which,  using  0  =  t/T,,  may  be  written: 

T, 

- .  0 

2  To 

B  —  Re  y  =  B  .  sin  (0  +  0 )  (47) 


41  is  obtained  by  addition  of  the  squared  equations 
37  and  39; 


equation  37 

42  is  the  quotient  of - 

equation  39 


We  see  that  the  amount  of  shut-down  has  influence 
on  the  size  of  R ,  but  not  on  the  size  of  0.  With  respect 
to  the  occurrence  of  the  movements  (see  equations  32  and 
33)  we  find  that  the  movement  is  a  damped  oscillation  with 
a  duration  of  8  secs  =  2»T,. 

Maximum  and  minimum  values  of  s  or  y  occur 
dz 

when  —  =0. 
dt 

This  is  the  case  when  sin  (7  —  0  —  t/T,)  =  o.  See 
Equation  35. 
t 

■ —  =  7  —  0=7  —  0  +  n  =  7  -  0  +  2*=  etc.  (43) 

T, 


The  former  gives  in  the  polar  system  of  co-ordinates 
(Fig.  3)  (0  in  radians)  with  0  as  the  vector  radius  a 
logarithmic  spiral  with  the  slope  “  =  —  TJ2T0. 

The  second  equation  may  be  demonstrated  in  the 
rectangular  system  of  co-ordinates,  where  y  and  there¬ 
fore  z  is  a  function  of  0,  if  we  consider  the  axis  of 
abscissae  in  the  continuation  of  this  ray  of  the  polar 
system,  which  is  turned  from  the  initial  point  with  the 
angle  0,  and  if  we  plot  from  an  initial  point  in  that 
system  the  values  of  0  projecting  the  points  of  the  spiral 
on  the  ordinates  which  apply  to  the  same  values  of  0. 
t 

As  0  =  — ,  the  time  may  be  measured  as  abscissae  show- 
T, 

ing  the  water  surface  fluctuations  with  respect  to  time. 

The  equation  35,  which  determines  5,  has  the  same 
form  as  equation  32.  It  follows  that  for  the  graphical 
demonstration  of  s,  the  same  method  may  be  applied  as  v 
becomes  o  when  </7’,  =7  —  0  and  we  see  that  the  axis 
of  the  abscissae  for  the  rectangular  system  of  co-ordinates 
which  demonstrates  s,  is  in  the  continuation  of  that  ray 
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of  the  polar  system  which  corresponds  to  the  maximum 
and  minimum  values  of  z  and  also  of  y.  By  projection 
of  the  points  of  #he  spiral  to  the  corresponding  ordinates, 
we  get  the  z  curve,  whose  ordinates,  as  may  readily  be 
seen,  cannot  be  measured  by  the  same  scale  as  those  of 
the  z  curve.  We  have  to  take:  If  (for  the  z  curve)  an 
ordinate  of  one  inch  equals  m  feet,  then  for  the  s  curve 

m 

an  ordinate  of  one  inch  =  —  feet  per  second.  The 

T 

demonstration  of  the  acceleration  is  similar. 


Vi  =  6.63  feet  per  second  and  h ,  =  9.6  feet,  are 

determined. 

The  surge  tank  section  may  be  constant  and  A  = 
5,380  square  feet.  According  to  the  theory,  we  assume 
hi  9.6 

that  h  =  n.v,  n  —  —  =  1.45  sec.  =  - .  Therefore, 

v,  6.63 

except  for  the  conditions  v  =  0  and  v  =  vt,  the  friction 
is  greater  than  it  would  be  according  to  the  friction 


formula  commonly  used  h  - 


L.a 


If  we  take 


also  applies.  The  dimensions  are  the  following : 

L  =  9,050  feet  p  —  32.8  feet 

a  =80  square  feet  0  =  530  cubic  feet  per  second 

For  a  flow  of  530  cubic  feet  per  second,  according  to  the 
usual  method  of  computation,  a  friction  coefficient  K  = 

L.a 

101,000  in  the  formula  h  =  - .  v3  and,  therefore, 

K.p 


R 2  sin  20  =  92  tgP  = 
R3  cos!0  =  1 19 


9.6 


=  —  .881  B  —  —  410  24' 


R2  =  211 


10.9 


2  To 


arc  B  = 


722 


tgy  = 


389 


147 


=  2.65  7  =  69°  18' 
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R  =  14.54  feet  arc  7  =  1.21 

t 


389  t 

z  =  14.54  e  sin  ( - .  722)* 

147 

t 

389  t 

s  =  .106  e  sin  (1.932 - ). 

147 

The  values  (in  seconds)  for  the  first  occurrence  of  the 
maximum  and  minimum  are 

t,  max  =  (7  —  0)  T,  =  [1.21  —  ( — .720)]  147  =  284  sec. 
f,  min  =  (7  —  0  +  *)  T,  =  (1.93  +  3.14)  147  =  746  sec. 

z  max  I.  =  6.55  feet,  z  min  II.  =  —  2.00  feet, 

8  --  2*Ti  =  924  sec.  =  15  min.  24  sec. 

The  graphical  demonstration  using  these  values 
shows  very  clearly  the  water  level  fluctuation  in  the  surge 
tank  in  function  of  the  time  (Fig.  3  and  Fig.  4). 

(4)  Work  Done. — For  the  case  of  a  sudden  complete 
shut-down,  it  is  interesting  to  strike  a  balance  of  the 
work  done.  The  basis  for  such  a  balance  may  be  secured 
by  a  variation  of  the  principal  equation  16,  considering 
also  the  continuous  equation  17.  If  we  multiply  the 
equation  16  by  (w  .a  .v  dt)h  we  get 

w  a  L 

- v  dv  +  (wav.dt)  z  +  (wavdt)h  =  o 

g 

w.a.L 

-  =  M  =  mass  of  the 

g 

main  conduit  volume  and  for  e  —  o,  w.a.v.dt  becomes 
=  w  A  s  dt  =  w.Adz  =  dG,  i.e.,  the  increase  of  weight 
of  the  contents  of  the  surge  tank  during  the  time  dt.  Also 
(w.a.dtv)  h  =  A1  =  to  the  amount  of  friction  work  in 
the  main  conduit  during  the  time  dt.  The  above- 
mentioned  equation  becomes 

dv 2 

M  -  +  zdG  +  dAl  —  o 

2 

and  the  integration  between  the  limits  which  refer  to  the 
initial  condition  and  to  the  highest  elevation  for  z  max., 
considering  that  at  the  beginning  v  =  v,  and  at  the  end 
v  —  o\ 

Mv*  C  +  2  max 

-  /  z  d  G  +  A'  =  o  -  -  (49) 

2  J  —h, 

We  may  place 

+  z  max 
z d  G  =  G.z, 

—  h, 

if  we  consider  G  as  the  weight  of  the  water  that  has 
entered  the  surge  tank  during  that  period  and  if  z,  is 
the  distance  of  the  centre  of  gravity  of  G  above  static 


*The  proportion  t/T ,  is  in  radians,  but  must  be 
changed  into  the  usual  measurements  as  soon  as  we  use 
the  trigonometric  functions. 


level  n  —  n.  But  we  may  as  well  say  Zi  =  yi  —  hi  (50) 
\y,  —  distance  of  the  centre  of  gravity  of  G  from  the  initial 
water  level).  It  follows  then,  that 

MV 

-  +  G  hi  —  Gyi  A1  -  (51) 

2 

* 

This  is  the  equation  for  the  work  done,  which  is  correct 
also  when  the  cross-section  A  is  variable.  As  potential 
energy,  we  must  introduce:  (1)  The  kinetic  force  of  the 
volume  of  the  main  conduit,  and  (2)  The  potential  energy 
of  the  weight  G  due  to  the  first  difference  of  level  hi  be¬ 
tween  the  elevation  in  front  of  the  forebay,  (i.e.,  the  static 
level  n  — n)  and  the  initial  level. 

As  used  work  we  must  introduce : 

(1)  The  lifting  work  due  to  the  lift  of  the  weight  G 
through  the  height  y,. 

(2)  The  friction  work. 

We  may  strike  the  balance  in  the  following  manner : 


Useful 

Used 

work  in 

work 

foot-tons. 

s 

15800 

Kinetic  Energy  of 
the  Conduit  Contents 

62  5  x  80  x  9050  6.63  = 

32.2  2 

26600 

Potential  Energy 
of  G  =  2770  Tons 
for  hy  =  9-6  Feet 

5380  (9.6  +  6.55)  x  62.5  x  9.6 

Lifting  Work  for 

2770  Tons  for  height 

16  15 

of  8.07  Feet 

22400 

2 

Friction  work 

20000 

42,400  —  22,400 

42400 

42400 

The  total  friction  work  of  20,000  foot-tons  corre¬ 
sponds  to  an  average  friction  head  corresponding  to  the 
conveyed  weight  of  2,770  tons. 

20,000 

hi  (average)  =  -  =  7.22  feet  =  .752.71,  (52) 

2,770 

We  might  use  this  result  to  determine  the  necessary 
surge  tank  volume  in  a  similar  case,  if  in  the  equation  51 
we  write:  A1  is  similar  to: 

h  (average)  .  G  which  must  be  equivalent  to  .752  /i,  G 

A  supplementary  formula  to  determine  the  value 
h  (average)  follows  in  the  final  chapter. 

Case  B. — Sudden  Opening. 

Before  the  beginning  of  outflow  through  the  pen¬ 
stock,  the  water  in  the  main  conduit  and  in  the  surge 
tank  is  quiescent.  The  water  surface  in  the  surge  tank 
is  at  the  level  n  —  n  and  z  =  zero.  As  soon  as  the  open¬ 
ing  occurs,  which  is  considered  as  taking  place  instan¬ 
taneously,  the  outflow  begins.  However,  the  full  amount 
of  the  outflow  does  not  occur  instantaneously  because  ol 
the  inertia  of  the  mass  ;  also  during  the  time  following 
the  opening,  q  will  not  be  constant  even  with  a  constant 
opening,  as  with  the  variation  of  the  water  surface  eleva¬ 
tion  the  slope  due  to  the  outflow  velocity  also  varies.  As 
in  case  (A),  where  we  assumed  that  all  these  influences 
and  variations  are  a  very  small  amount,  we  consider  q  as 
constant  from  the  instant  that  the  opening  occurs. 

The  following  will  now  happen  at  this  juncture:  Im¬ 
mediately  after  the  opening  the  water  surface  in  the  surge 
tank  goes  down.  The  head,  due  to  this  difference  between 
the  water  level  in  the  main  conduit  and  that  in  the  surge 


IO 


SURGE  TANK 


P  R  O  B  L  EMS. 


tank,  creates  a  movement  in  the  main  conduit,  but  on 
account  of  the  inertia  of  the  mass  and  the  increasing 
friction,  the  velocity  does  not  increase  in  the  same  pro¬ 
portion  as  the  drop  of  the  surface  in  the  surge. tank  below 
the  static  level  n  —  n.  On  this  account  the  water  surface 
in  the  surge  tank  goes  below  the  level  which  corresponds 
to  the  normal  water  level  for  a  constant  flow.  So  that  we 
have  first  a  dropping  down  to  a  lowest  level  whereas 
afterwards  there  occurs  a  damped  oscillation,  until  finally, 
under  the  influence  of  all  resistances,  a  normal  balanced 
condition  is  reached,  where  the  level  in  the  surge  tank  is 
as  much  under  the  level  n — n  as  the  amount  of  head  re¬ 
quired  to  create  the  necessary  slope. 

The  initial  equation  23  has  in  this  case  (if  we  say 
q  =  €0,  for  purposes  of  comparison)  the  same  form  as  in 
case  (.4)  and  therefore  all  other  equations  (growing  out  of 
equation  23)  also  apply.  The  determination  of  the  in¬ 
tegration  constant  will,  however,  give  other  values  for 
R  and  P  since  the  initial  condition  is  different.  We  have 

0. 

at  the  time  2  —  zero,  z  =  zero,  and  s  =  s„  —  —  « —  = 

A 


mentioned  investigations)  it  follows  that  the  initial  radius 
vector  in  the  polar  co-ordinate  system  for  the  case  ( B )  is 
turned  through  1800  from  case  (A).  The  duration  of  the 
period  is  the  same  as  in  case  (A).  The  graphical  demon¬ 
stration  may  be  begun  also  from  the  same  spiral  with  the 
slope  tg  «  =  — ■  T 1/2 To.  The  analytical  results  (assuming 
a  full  opening  from  zero  to  530  cubic  feet  per  second, 
that  is,  6  =  1,  n  —  1.45  sec.,  assuming  the  same  di¬ 
mensions)  are  as  follows : 

T  —  137.5  sec.  To  =  194.5  sec-  Ti  =  147  sec. 

R  —  14.54  feet  0  =  180  +  ( —  41 0  24')  =  138°  36' 

arc  P  =  +  2.420 
—  2/389 

z  =  —  9.6  +  14.54  e  sin  (2.42  +  f/147) 

y  =  69°  18'  arc  y  —  1.210 
—  P38Q 

5  =  .106  e  sin  ( —  1. 21  —  2/147) 

—  2/389 

=  —  .106  e  sin  (1.932  —  2/i47) 

Therefrom  we  determine  the  time  for  the  first  mini¬ 
mum  of  z 


Fig.  4. 


—  «.Ci.  For  the  determination  of  the  constants  R  and  P 
we  get  the  equations  : 

ht 

R .  sin  P  =  +  *hi ;  R  cos  P  =  —  e(C, - )  T 1  (53) 

2  To 

The  same  procedure  as  in  case  (A),  referring  to  the  signs 
of  sin  P  and  cosP,  shows  in  this  case  that  P  lies  in  the 

A  2  L 

second  quadrant,  if  —  <  - ,  and  in  the  first  quadrant 

a  n2g 

4  L  A  2  L 

when  -  >  —  > - . 

n2g  a  n'g 

We  get,  further, 

T.To  1 

R  =  *h, -  tg  P  = - -  (54) 

T2  To/T»  —  }  Ti/T„ 

If  we  open  the  penstock  so  that  the  outflow  has  the 
full  value  q  =  Q,  (e  =  1)  then  R  assumes  the  same  value 
as  in  case  (A)  for  a  complete  shut-down.  P  has  in  both 
cases  the  same  form,  therefore,  (referring  to  the  above- 


2.  min  =  1.932  •  147  =  284  sec. 
and  for  the  first  maximum : 

2,  max  =  (1.932  +  3.142)  147  =  746  sec. 

These  are  the  same  values  as  were  reached  before  for 
the  occurrences  of  the  first  maximum  and  minimum. 
Finally  follows 

z1  min  —  —  9.6  —  6.55  =  —  16. 15  feet ; 
z'max  =  —  7.6  feet. 

As  z'  max  is  negative,  we  see  that  with  the  as¬ 
sumed  dimensions  and  for  full  opening  to  530  cubic  feet- 
seconds,  the  water  surface  does  not  rise  any  more  above 
the  level  n  —  n  after  the  first  dip.  (Fig.  4.) 

The  velocity  in  the  main  conduit  at  the  time  of  the 
greatest  drop  becomes  Vi  =  QJa  but  the  velocity  increases 
continuously  thereafter,  as  in  the  following  period  of  rise 
of  the  water  surface,  besides  the  outflowing  quantity  Q, 
also  that  water  from  the  main  conduit  which  is  necessary 
for  the  refilling  of  the  surge  tank  has  to  flow  in.  And 
the  increase  lasts  until  v  becomes  a  maximum.  Because 
v.a  =  s  A  +  0 1  the  velocity  v  becomes  a  maximum  if 
x  is  a  maximum  and  we  get : 
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*t2 


R  2  Ta 

s  max  =  — e  .sin  y  =  .03  feet  per  second. 

T 

(55) 

A 

v  max  —  —  .  s  max  +  Vi  —  8.65  feet  per  second. 

a 

The  conclusion  of  practical  value  is  that  the  greatest 
drop  below  the  static  level  n  —  n,  resulting  from  a  sudden 
full  opening,  reaches  the  same  value  as  the  greatest  rise 
from  the  initial  level  for  a  sudden  full  shut-down. 

Case  C:  Variable  Outflow. 

In  the  following,  we  consider  first  the  influence  of  a 
gradual  shut-down  and  opening  of  the  penstock,  and 
secondly  the  case  when  the  outflow  is  variable  with  re¬ 
spect  to  time.  The  movement  during  the  period  of  a 
gradual  shut-down  will  be  quite  different  from  that  after 
the  shut-down.  Therefore,  we  must  consider  both  cases 
separately. 

(1)  Movement  During  the  Shut-down. — We  assume 
that  the  closing  occurs  so  that  the  outflow  decreases  in 
direct  proportion  to  the  time  which  is  expressed  by  the 
formula  : 

t 

q  —  Q'  (1 - )  (56) ;  t  =  time  for  complete  shut-down 

T 

Therefore, 

t  dc  Ci 

C  =  Ci  ( 1 - )  —  = - 

T  dt  T 

Similar  to  equation  (23)  we  get  the  special  principal  equa¬ 

tion  (57)  in  the  following  form  : 
d2z  1  dz  z  Ci  t  Ci 

-  +  —  .  —  +  —  - (1 - ) - =0  (57) 

dt 2  To  dt  T2  To  t  t 

The  general  integral  of  this  equation  may  be  obtained 
by  means  of  the  theory  of  linear  differential  equations  of 
the  second  order. 

t 


2  To  t 

The  general  member  is  Zi  —  R  e  sin  (0  +  — ) 

T, 

The  particular  integral  z2  =  hi  +  b2.t  when  3,  and 
b2  may  be  determined  by  inserting  the  value  of  z2,  also 
dz  d2z 

—  =  b2  and  -  =  o  in  equation  57. 

dt  dt2 

Therefore, 

b  2  b  1  b2t  Ci  cd  Ci 

—  +  —  +  -  + - =0  -  (58) 

To  T2  T2  To  To  t  t 

As  the  conditions  of  this  equation  must  be  fulfilled  for  all 
values  of  t  the  two  equations  follow: 

b  2  b  1  Ci  Ci  b2  Ci 

-  +  —  + - =  0;  — - =  0  (59) 

To  T2  To  t  T2  Tot 

ClT2  To  T2  hi 

As  hi  =  -  b,  =  hi  [ —  (1 - )  —  1]  b2  —  — 

To  t  To2  t 

The  general  integral  is  obtained  by  the  form  z  =  z,  +  z2. 


Therefore, 

3  =  R.e 


- t/2To 

sin  (0  +  t/T  1)  +  hi  +  b2t 


And  with  differentiation  with  respect  to  t  (tg  y  = 
R  — t/2T0 

s  —  —  e  sin  (7  —  0  —  t/T  1)  +  b2 

T 


(60) 

2  To 


Ti 

(61) 


In  order  to  obtain  the  constant  for  the  integration, 
consider  that  for  t  =  o,  z  =  Zo  =  —  hi,  as  in  case  (A),  but 
that  s  =  So  =  o,  because  the  shut-down  occurs  gradually. 
Therefore : 

T„  T 

R  sin0  =  —  hi  —  [1  - - ] 

T  To2 


Ti  T2 

R  cos  0  =  —  hi - [3 - ] 

2T  To2 

T2 

1 - 

T„  Ti  T„2  T„ 

R  —  hi  .  —  .  —  (62)  tgP  =  2  -  .  —  (63) 

T  T  T2  Ti 

3 - 

To2 

For  the  graphical  demonstration  of  these  functions 
we  may  separate  them  into  three  equations.  For  instance, 
the  s  —  function  into 

- tl2To<P 

r  —  Re  ;  3i  =  r  sin  (0  +  ?>);  =  hi  +  {b2Ti)  <t>.  (64) 


The  value  of  r  determines  again  a  logarithmic  spiral 


Ti 

with  the  slope  tg  a  = - with  the  initial  vectors 

2  To 

r0  and  0. 

3,  is  obtained  in  the  rectangular  co-ordinate  system 
by  projection  from  the  polar  system  as  in  the  former 
cases.  z2  in  the  same  rectangular  co-ordinate  system  is 
a  straight  line  which  intersects  the  ordinate  axis  at  the 
distance  ht  from  the  initial  point  and  whose  inclination 
to  the  axis  of  abscissae  is  fixed  by  the  direction  constant 
b2.Ti.  The  algebraic  sum  of  zi  +  z2  gives  z.  (See  Fig.  5.) 
But  only  that  part  of  the  curve  which  lies  between  the 
values  of  the  abscissae  o  and  t  is  of  practical  importance 
because  at  the  time  t  the  complete  shut-down  has  already 
occurred  (as  we  assumed). 

(2)  Movement  After  Completed  Shut-down. — After 
the  shut-down  has  occurred,  the  movement  of  the  water 
surface  continues  as  described  in  case  (A).  Therefore, 
from  this  moment  the  following  equations  are  effective : 


—  R2T0 

z  =  R,e  sin  (0,  +  t/T,)  -  -  (65) 

R,  — tj  2T0 

s  =  — e  sin  (7  —  0.  —  f/T,)  (66)  (see  equations 

T  32  and  35). 

2  To 

Where  tg  7  =  - 

.  T, 

The  constants  /?.  and  0.  have  to  he  determined  from 
the  value  of  3  and  s  computed  for  t  =  t.  In  order  to 
obtain  an  easier  computation,  we  measure  the  time  for 
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the  second  period  from  the  moment  of  the  completed  shut¬ 
down,  so  that  the  following  equations  for  the  determina¬ 
tion  of  the  constants  R i  and  /3,  are  in  effect : 

R i  sin  0,  =  zT 

zT 

Ri  COS/3,  —  (5t  4.  - ) 

2  To 

where  zT  and  sT  are  limiting  values  of  the  first  period. 


(3)  Practical  Example. — The  computation  in  con¬ 
nection  with  the  preceding  example  follows  for  a  shut¬ 
down  in  10,  100  and  200  seconds,  with  a  discharge  of  530 
cubic  feet  per  second. 

From  the  graphical  demonstration  (Fig.  4)  we  see 
that  the  period  of  shut-down  up  to  100  seconds  influences 
the  maximum  elevation  a  very  small  amount  indeed.  The 
explanation  for  this  is  that  during  this  period  the  velocity 


The  demonstration  of  the  movement  may  be  developed 
according  to  the  same  method  as  in  case  (d).  For  the  com¬ 
putation  of  the  value  of  the  most  practical  importance 
(the  maximum  elevation),  we  get 

T, 

- —  (7  —  /3.) 

2  To 

sin  7  (67) 


of  flow  in  the  main  conduit  decreases  very  little.  For  the 
determination  of  the  dimensions  of  the  surge  tank,  we 
have,  therefore,  to  consider  the  results  of  the  limiting 
case,  that  is,  the  sudden  shut-down.  The  computation  for 
that  condition  is  a  much  simpler  one,  (case  .4).  For  a 
gradual  opening,  the  same  relations  exist.  The  lowering 
of  the  level  n — -n  occurs  to  the  same  extent  and  for  the 
same  duration  as  the  rise  above  the  initial  level  in  the 
case  already  computed. 


1 


0  max  =  Ri  e 


CASE  C:  VARIABLE  OUTFLOW. 


*3 


■# 


i 


With  T  =  137.5  seconds,  To  —  194.5  seconds,  7\  =  147  seconds,  hi  —  9.6  feet,  we  get  the  following 
table : 


for  t 

= 

10  sec. 

100  sec. 

200  sec. 

To 

T3 

bt  =  hi  [ —  (1- 

+  83.3  feet 

— .264  feet 

—  4.93  feet 

T 

To 2 

hi 

feet 

feet 

feet 

b2  =  — 

= 

.96 

.096  - 

.048 - 

T 

sec. 

sec. 

sec. 

To 

Ti 

R  =  hi  .  —  . 

—  = 

200  feet 

20  feet 

10  feet 

T 

T 

T2 

2  (i  - 

-—)  To 

To3 

tg&  - 

-  = 

•53° 

•530 

•530 

T2 

(3  — 

— )  Ti 

To3 

0  =  207°  54' 

arc  0  =  3.628 

7  ==  69°  1 8' ;  arc  7=1.21 

zT  = 

s  T  = 

Ri  sin  0,  = 
R»  cos  0i  - 
R.  = 

0,  negative  = 
arc  0i  = 
3  max  = 


— 9. 10  feet 
feet 

.097 - 

sec. 

— 9.  xo  feet 
+  1 1.05  feet 
14.30  feet 
—39°  30' 
— .6894 
+  6.53  feet 


— 4.85  feet 
feet 

.091  - 

sec. 

— 4.85  feet 
+  1 1.6  feet 
12.55  feet 

- 22°  44' 

—.3968 
+  6.40  feet 


— 1.05  feet 
feet 
•0735  — 
sec. 

— 1.05  feet 
+  10.4  feet 
10.45  feet 
-5°  46' 
— .  1006 
+  5.97  feet 


Of  special  interest  is  the  case  of  an  outflow,  which 
is  variable  in  the  sense  that  the  outflow  increases  con¬ 
siderably  during  a  certain  time  and  decreases  afterwards 
to  the  same  amount  as  before  or  to  some  other  amount, 
(for  instance,  in  a  plant  for  railway  operation). 

x.  Analytical  Investigation. — We  assume  now,  that 
under  circumstances  similar  to  those  above  mentioned, 
the  following  law  of  outflow  is  effective: 

q  =  *-Q  1  (t  +  /  sin  t/ t) 
so  that  for  the  time 


and  further, 

1  t  1  t  1 

e./.Ci  [ —  sin  —  +  —  cos  —  ]  = 

T„  T  T  T 


=  ./.c  \/  —  +  —  si 

V  To3  t3 


—  +  —  sin  (<P  +  — ) 


with  tg  <P  =  —  then  the  equation  69  transposes  to 

T 


(7i) 


5T.T 

t  =  0  :  t  =  — ;  t  =  q  becomes  resp.  =  *Q,  ; 

2 

*Q,  (1  +  /);  -  -  (68) 

after  the  time  t  =  * . T  the  outflow  may  be  constant  again. 

e  and  /  are  natural  numbers.  /  is  the  proportion  of 
the  maximum  increase  of  the  outflow  to  the  normal  out¬ 
flow.  So  that 

q  t  dc  Ci  t 

c  =  —  =  «ci  (1  4-  /  sin  — )  —  =  «./.  —  cos  — 

A  t  dt  T  T 

Equation  23  may  then  be  written : 


d3y 

dt3 


+ 


1  dy  y 

- +  -  + 

T„  dt  T3 


1  t 

+  —  sin  ( <P  H - ) 

T*  T 


o  (72) 


and  it  then  follows  again  from  the  theory  of  differential 
equations  of  the  second  order  that  y  =  yi  +  y=,  that  is, 

-  t/2To 

y,  =  R.e  sin  (0  +  t/T,)  general  integral  (73) 


d3Z  I  dz  Z  «.C  1  e. /.Cl  t  e./.Ci  t 

- 1 - . - 1 - 1 - 1 - sin - 1 - cos  —  =  0  (69) 

dt3  To  dt  T3  To  To  T  T  t 

If  we  introduce 

« Ci  T3  dz  dy  d3z  d3y 

2  =  y - =  y  —  fh, ;  —  =  — ; - = -  (70) 

To  dt  dt  dt3  dt3 


t 

y2  =  b  sin  (i>  +  — )  particular  integral  (73) 

T 

The  latter  may  also  be  written 

t  t 

y2  =  b.  sin  ^  cos  —  +  b  cos  $  sin  —  (74) 

T  T 
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The  integration  constants  R  and  0  may  be  obtained 
from  the  initial  phase  and  the  constants  b  and  *p  by  form¬ 
ing  the  equations 

dy2  b  t  b  t  b  t 

—  =  —  cos  (V'  +  — )  =  —  cos  'P  cos - sin  41  sin  — 

dt  T  T  T  T  T  T 

d2y2  b  t  b  t  b  t 

—  = - sin  ('P  H - )  = - sin  'P  cos  — - cos  'P  sin  — 

dt 3  T2  T  T2  T  T2  T 

We  then  introduce  these  values  in  the  differential  equation 
and  combine  the  members  containing  cos  t/r  and  sin  t/ t 

b  b  b 

[ - sin  +  - cos  'P  +  —  sin  'P  + 

t2  t  T0  T2 


,.c  J-  . 

>  7’.2  T2 


t 


—  +  —  .  sin  <p]  cos  —  + 
1\2  t2  t 


b  b  b 

[ — -  —  cos  'P - sin  t  +  —  cos  ^  + 

T2  T  To  T2 


I  t 

+  —  .  cos  <P]  sin  —  —  o 

T2  T 


and  as  these  equations  hold  good  for  all  values  of  t,  the 
terms  in  both  parentheses  must  become  o.  One  thus 
obtains  two  equations  with  the  unknowns  b  and  'P  and 

To  Ci  T2 

considering  that  tg<p  =  —  and  hi  =  - ,  it  follows 

T  To 


b 


o/.fc, 


tg+  = 


rp2  rp2 

I  +  -  +  —  2] 

T2  T2  T„2 

To  T2  T2 

--•[-  + - 1] 

T  T2  To2 


(75) 


For  the  determination  of  the  integration  constants 
R  and  0,  it  must  be  considered  that  for 


t  =  o  ;  z  =  —  ;  y  =  o  ;  s  =  o. 

We  obtain  for  this  case  the  equations 
R .  sin  0  ==  —  b  .sin  'P 

Ti  Ti  (76) 

R  ■  cos  0  =  —  b  [ —  cos  *P  +  - sin  'P] 

T  2  To 

During  the  variation  of  the  outflow  the  following  equa¬ 
tions  are  effective : 


t 


■\ 


2  T„  t  t 

z  —  ^.hi  +  Re  sin  (0  H - )  +  b  sin  ('/'-) - 

Ti  T 

t 


with 
2  To 

tgy= - 

Ti 


R  2  To  t  b  t 

s  —  — e  sin  (7 — 0 - ) -) - cos  ($1 - ) 

T  T,  T  T 


> P  amounts  generally  to  about  —  which  relation  simplifies 

2 

the  computation  of  R  and  0. 

The  variation  ceases  (as  we  assumed)  after  the  time 
t  =  n-.T.  The  values  which  correspond  to  the  elevation 
of  the  water  surface  and  velocity  at  that  time  are  obtained 
by  the  formulae : 

T 


2  T„  T 

=  —  *hi  +  R  e  sin  (0  +  «■  — )  —  b .  sin  ^ 

Ti 

T 


R  2  To  t  b  (77) 

sT  =  —  e  sin  (7  —  0  —  *  — ) - cos  'P 

T  T,  t 

For  further  investigation,  the  formulae  of  case  (4) 
may  be  applied.  In  order  to  simplify  matters,  in  the  de¬ 
termination  of  the  further  movement,  we  measure  the  time 
from  the  instant  of  the  beginning  of  the  constant  outflow 
and  logically  we  must  use  the  limiting  values  of  the  pre¬ 
ceding  phase  for  the  determination  of  the  integration  con¬ 
stants  Ri  and  0i. 

If  we  do  not  hinder  the  variation  of  the  outflow  but 

t 

maintain  the  law  *Qi  (1  +  /  sin  — ),  we  see  that  the  move- 

T 

ment  of  the  water  surface  in  the  surge  tank  takes  the 
form  of  a  forced  oscillation  ;  where  the  influence  of  the 
first  member  decreases  with  the  increase  of  t  and  this  the 

Ti 

quicker  the  larger  the  value  of  —  becomes  in  the  member 

2  To 

t  Ti  t 

2  To  2  To  T, 

e  =  e 

The  movement  of  the  level  of  the  water  surface  is 
merely  that  of  a  harmonic  oscillation.  In  such  cases  it  is 
well  known  that  the  phenomenon  of  resonance  may  occur, 
if  the  period  of  the  actuating  influence  has  the  same  dura¬ 
tion  as  the  swinging  bodies’  own  period,  that  is  to  say, 
if,  in  the  case  mentioned  T  —  T.  The  value  of  the  ampli¬ 
tude  of  the  forced  oscillation  is  then 


We  obtain  the  angle  'P  from  tgi>  =  —  —  = - (79) 

To  tg<P 


To 

because  —  =  tg<p  and  t  =  T. 

T 

According  to  the  theory  of  the  phenomena  of  reson¬ 
ance,  the  amplitude  of  the  forced  oscillation  becomes  of 
infinite  value  when  the  damping  force  is  infinitely  small. 
The  latter  would  be  the  case,  using  our  terms,  when  To  = 
infinite  with  which  value  the  amplitude  of  the  above- 
mentioned  equation  also  becomes  infinite.  The  phase 
difference  between  force  and  movement  has  thereby  to  be 
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i5 


equal  to  — ,  which  also  finds  its  expression  in  this 
2 

equation  (79). 

2.  Graphical  Demonstration. — The  graphical  de¬ 
monstration  of  the  movement  of  the  first  phase  is  not  dif¬ 
ficult.  The  demonstration  of  the  values 

t 


2  To  t 

Zi  =  R  e  sin  (ft  4-  — ) 

Ti 

may  be  carried  out  in  the  same  form  as  before  by  means 
of  the  projection  ©f  the  logarithmic  spiral.  The  demon¬ 
stration  of 


—  .  t  =  31.4  seconds  to  530  cubic  feet  per  second,  and 
2 

decreases  after  another  period  of  31.4  seconds  (a  total  of 
62.8  seconds)  to  265  cubic  feet  per  second.  Therefore,  with 
T=  137.5  sec.  T0=  194.5  sec.  Ti=i47sec.  ^  =  9.6  ft. 

To  T  T 

we  get  —  =  9-725:  —  =  6.875;  —  =  -707 
T  T  To 


To2  T2  T2 

!+-[-  +  -  -I]* 

T2  T2  To2 

b  =  —  efh ,  -  -  =  —  1.008  feet 

rp2  J2  rp2 

1  +  —  [—  + - 2] 

T2  T2  T„2 


is  obtained  by  projection  of  a  circle,  constructed  with  the 
radius  b  in  the  polar  system,  where  the  initial  vector  is 
inclined  against  the  axis  of  abscissae  with  an  angle 

In  this  case  a  rotation  of  the  initial  vector  of  z,  by 
an  amount  of  2  *  corresponds  to  a  rotation  of  the  initial 

Ti 

vector  of  z2  by  an  amount  equal  to  —  21.  (See  Fig.  6.) 

T 

3.  Analytical  Example. — We  use  the  same  conduit 
and  surge  tank  dimensions  and  assume  the  same  friction 
conditions.  As  before,  we  take  «  =  0.5;  with  f  —  1 
and  for  t  =  20  seconds  we  have  an  initial  flow  of  265 
cubic  feet  per  second  which  increases  during  the  time 


R  =  b 


✓ 


T,2 


1  + 


(2  T„)2 


=  b  .7 


•V-. 

V  T,2 


(2  To)2 


T,  147 

b  .  —  =  —  1.008  . - 

T  137-5 


R  —  —  1.075  feet>  therefore 
t 


25  =  —  4.8  —  i  ,075 .  e 


389  t  t 

sin  [69°  18'  -I - ]+  1.008  cos  — 

147  20 
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t 


389  t  t 

v  =  +  .  00782  e  sin  —  —  .0503  sin  — 

147  20 

If  at  the  time  t  =  i.t  the  fluctuation  of  the  outflow 
ceases  it  follows  that 


/  =  1  the  forced  oscillation  becomes  nearly  a  constant 
harmonic  oscillation  of  the  amplitude  one  foot  measured 
from  2  =  —  4.8  feet  and  for  a  duration  of  period  of 
137.6  seconds. 

The  case  of  resonance  might  occur  if  T  =  T  =  137.5 
seconds.  This  would  be  the  case  for  an  amplitude  of  the 
forced  oscillation  (measured  from  z  =  —  4.8  feet)  equal  to 


20 

389  20 


ZT  — 

—4.8— 1.075  e 

sin  (69°  1 8' 

N 

II 

O 

i 

.88  ft. 

T47 

20 

389 

20 

vT 

=  .00782  e 

sin  —  =  .00275  feet  per  second. 

*47 

For  the  movement 

:  that  follows, 

we  may  use  the 

equa- 

lions  (32)  and  (35) 

t 

2  To 

t 

z  = 

—  4 -  Rx  e 

sin  (0i  4- 

-) 

t 

Tx 

R,  2  To 

t 

2  To 

s  = 

—  e  sin  (7  —  0i - ) 

with  tg  7  = 

— 

T 

Tx 

Tx 

and 

R x  and  0i  may  be  determined  by 

the  conditions  t 

=  Ox 

z  = 

aTi  5  =  5 T.  We  obtain 

Rx  . 

sin  0i  =  —  9.6  4- 

2.88  4-  4.8  ~ 

--  —  1.92  feet 

Rx  . 

147 

cos  0i  =  .00275 

ON 

1 

.32 

1  feet 

389 

Rx  = 

=  1.95  feet ;  0,  lies 

in  the  third  quadrant,  i.e.,  negative. 

tg  0 

*  =  5-98975 

-  99°  3°' 

Arc  0i  =  — 

i-737 

7  = 

68°  i87  arc  7  = 

=  1,210 

t 

389 

[—  99°  3°'  + 

t 

2  - 

—  4-8  +  i-95  e 

.  sin 

- ] 

T47 

J  T„2  To 

b  j=  ±  e./.fri  4  /  - - b  i.  —  =  ±  11. 7  feet. 

V  t2  T 

Herein  the  duration  of  period  would  be  T.  21  =  137.5.2.1 
sec 

=  863  =  i4;  23" .  It  does  not  seem  impossible  for 

such  periods  to  occur  in  regular  operations  as  for  rail¬ 
ways,  so  that  this  investigation  should  be  useful. 

The  preceding  method,  introducing  a  periodical  func¬ 
tion  for  the  determination  of  a  variable  outflow  in  function 
of  the  time  may  be  enlarged,  introducing 

t  t 

q  —  €£?i[/o  +  /i  sin  (0,-1 - )  4-  /2  sin  ($>2+2  — )  4-  .  .  •]  (80) 

T  T 

where  we  find  the  values  for  /  and  <t>  corresponding  to  a 
given  variation  of  q  by  means  of  the  Fourier  series.  The 
integration  of  the  differential  equation  thus  obtained  in¬ 
volves  no  difficulty  whatever.  It  is  based  upon  the  same 
method  as  that  given  in  the  last  example.  The  values  of 
2  and  5  have  the  form 

t 


2  To  t  t 

z  =  20  4-  R  -e  sin  (0  4-  — )  4-  2  fen  sin  ('/'n  4-  n — ) 

Ti  T 

t  (8l) 

R  2  To  t  n  .  sn  t 

s  =  Jo  +  —  e  sin  (7  —  0  4 - )  4-  2 - cos  (^n  4-  n  — ) 

T  T,  T  T 

Naturally  the  computation  requires  great  care.  The  gra¬ 
phical  demonstration  is  obtained  by  superposition  of  the 
projections  obtained  from  the  polar  system  with  the 

Tx  t 
2  To  Ti 


389  t 

s  =  —  .0142  e  .  sin  [n°  i27  4 - ] 

147 

Since  at  the  time  t  —  o  the  velocity  5  is  positive,  a  maxi¬ 
mum  first  occurs.  The  time  from  the  beginning  of  that 
maximum  until  the  beginning  of  the  second  phase  is 
tx  =  147  (7  —  0i)  =  433  seconds. 

433 


389 

z,  min  =  —  4.8  4-  1.95  e  sin  7  =  —  4. 19  feet 

The  general  course  of  the  movement  may  be  seen  in 
Fig.  6. 

In  case  of  a  continuous  fluctuation  of  flow,  expressed 

t 

by  the  law  q  =  e0i  (1  +  f  sin  — )  with  *  =  .5  and 

T 


logarithmic  spirals  r  =  R  e  and  the  circles 

with  the  radii  k. 

Case  D:  Involving  Spillways. 

In  order  to  decrease  the  surge  in  the  surge  tank,  a 
spillway  may  be  introduced  in  the  surge  tank  or  in  the 
main  conduit.  The  problem  of  interest,  then,  is  to  deter¬ 
mine  at  what  elevation  and  with  what  capacity  such  a 
spillway  should  be  designed  in  order  that  the  surge  shall 
not  exceed  a  certain  elevation.  The  case  will  be  investi¬ 
gated  under  the  following  assumptions :  The  spillway  is 
attached  to  the  surge  tank  (Fig.  7).  The  crest  of  the 
spillway  is  at  the  distance  e'  from  the  static  level  n  —  n; 
e'  is  positive  if  the  crest  is  above  n  —  n  and  negative  if 
below.  The  width  of  the  spillway  is  b1.  The  outflow  of 
Q  cubic  feet  per  second  is  suddenly  interrupted. 

We  must  take  care  of  several  periods  of  movement. 
Following  a  complete  shut-down,  the  water  surface  in  the 
surge  tank  will  rise  to  the  height  of  the  spillway  crest,  as 
is  described  in  case  (A).  Next,  the  rise  will  continue,  but 
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water  will  flow  over  the  spillway  up  to  a  maximum  eleva¬ 
tion,  then  decrease  to  the  elevation  of  the  spillway  crest. 
At  this  moment,  the  outflow  ceases.  The  movement  then 
takes  the  form  of  the  case  (A),  and  retains  this  form  there¬ 
after,  provided  the  crest  of  the  spillway  is  not  reached 
again.  Otherwise,  the  same  cycle  is  repeated  or  the  con¬ 


dition  remains  one  of  constant  overflow.  Therefore,  for 
the  first  phase  the  following  equations  are  effective: 
t 


2  To  t 

z  =  R.e  sin  (0  +  — ) 

T, 

t 


quantity  equals  zero,  which  is  the  same  as  if  we  move  the 
elevation  of  the  spillway  crest  higher  up  by  the  same 
amount,  and  if  we  take  from  this  location  the  water 
quantity  proportional  to  the  height  of  this  new  elevation 
of  the  spillway  crest  corresponding  to  the  tangent,  then 
we  get  a  linear  equation  for  the  determination  of  the 
water  quantity,  whereas  the  differential  equation  for  the 
determination  of  the  elevations  also  becomes  linear. 

The  results  thus  obtained  from  these  equations  are 
naturally  only  approximate  ones,  as  the  overflowing  quan¬ 
tities  are  introduced  as  too  small.  The  computed  values 
of  z  exceed  the  actual  values.  For  practical  purposes, 
this  first  approximation  is  generally  sufficient,  but  we 
have  no  difficulties  using  the  results  of  the  first  approxi¬ 
mation  for  a  second  computation,  drawing  the  tangent 
on  that  point  of  the  overflow  curve  which  corresponds  to 
the  maximum  value  of  the  elevation  found  in  the  first 
computation  and  repeating  with  these  results  the  compu¬ 
tation  as  before.  We  may  use  this  second  approximation 
in  the  first  computation  if  instead  of  the  value  for  the 
point  in  the  curve,  we  take  a  somewhat  smaller  value, 
say,  uAsA,  when  u  =  0.7  to  0.8.  We  use  the  latter 
method  in  the  following: 

The  overflow  height,  which  gives  an  overflow  of 
u.A  .sA,  is  determined  by 

3  U  sA  A 
hn  =  ( - .  - ——)*/■ 

2  A*  b1  A  2g 


R  2  To  t  2  To 

s  —  —  e  sin  (7  —  0 - )  with  tgy  —  - 

T  Tt  T, 

The  integration  constants  R  and  0  are  obtained  for  the 
conditions  t  =  o,  z  =  — ft,,  s  =  c  ,.  The  final  value  of  z  is 
in  this  period  of  movement  e'.  From  this  we  determine 
the  time  tA ,  which  is  necessary  for  the  rise  up  to  the  level 
of  the  spillway  crest  and  by  means  of  the  second  equation 
we  determine  the  final  velocity  sA.  In  the  second  period 
of  movement  now  beginning,  we  have  overflow  on  the 
spillway,  that  is,  according  to  the  familiar  formula  for 
spillways  q  =  f  m  b'  h'  A  2  g  h'  where  h'  =  z  —  e'  repre¬ 
sents  the  overflow  height  and  we  get 

q  A  2  ghnr 

—  =  c  =  t  a*  ft' -  (82) 

A  A 

The  introduction  of  this  formula  and  its  derivation  in  the 
principal  equation  (23)  would  lead  to  a  differential  equa¬ 
tion  higher  than  the  first  degree,  the  integration  of  which 
might  be  accomplished  by  development  in  series.  For 
practical  purposes,  however,  an  easier  but  sufficiently 
exact  approximation  may  be  obtained  in  the  following 
manner : 

If  we  plot  the  different  values  of  h'  computed  for  the 
width  b'  on  a  rectangular  co-ordinate  system,  the  values 
h'  as  abscissae,  and  the  values  of  q  as  ordinates  (see  Fig. 
8),  we  obtain  a  parabolic  curve,  passing  through  the 
origin.  The  velocity  sA  (which  is  the  final  value  of  the 
velocity  of  the  surface  at  the  end  of  the  first  period  of 
movement)  together  with  A  determines  a  quantity  of 
water  sA.A,  which  is  certainly  larger  than  the  maximum 
value  of  the  overflowing  quantity  during  the  second 
period.  If  we  now  draw  tangents  from  the  point  of  the 
quantity  curve,  which  corresponds  to  Ase;  and  assume  as 
a  preliminary  approximation  that  within  the  heights  of 
the  overflow,  which  corresponds  to  the  intersection  of  the 
tangent  with  the  axis  of  the  abscissae,  the  overflowing 


The  proportional  factor  k  for  the  linear  variation  of  q  is 
obtained  by  differentiation  of  q  with  respect  to  k. 
Therefore : 
dq 

k=  —  =  (as  h*  =  ha)  =  n  b'  A 2 g  ftu  =  ^3/2  /A  b'2  2 gu  sA A 
dhl 


k  has  the  dimension  and  the  value  of  the 

abscissa  e,,  which  is  the  difference  between  the  true  ele¬ 
vation  of  the  spillway  crest  and  that  obtained  by  ap¬ 


proximation,  is 

u.sA.A  f  19  u  sA  A 

e.  =  ft. - e,  =  J7  —  ( - -yr  (84) 

ft  T  12  m  b'  A 2 g 

The  values  are  easiest  obtained  graphically  from  the 
curve  of  the  overflow  quantities. 

Therefore,  with  e,  the  height  of  the  ideal  spillway 
crest  above  the  static  level  n  —  n,  (that  is  E  =  e1  +  e>) 
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once  determined,  the  computation  of  the  first  period  of 
movement  must  be  extended  to  the  elevation  E. 

We  obtain  from  the  previously  mentioned  formulae 
Ze  =  E  and  se.  These  are  initial  values  for  the  second 
phase,  from  which  beginning  we  measure  the  time  anew. 


by  which  formula  the  investigation  mentioned  may  be 
carried  out  and  the  corresponding  form  of  the  general 
integral  may  be  used. 

The  integration  constants  must  be  determined  with 
the  initial  values  y 


q  k  dc  k  dz 

c  =  —  =  —  (z  —  E)  and  therefore  —  = - 

A  A  dt  A  dt 

and  the  equation  23  becomes 

d*z  1  k  dz  1  k  k 

-  +  ( —  +  — )  —  +  ( —  +  - )  2  —  E . -  =  o 

dt 1  T„  A  dt  T*  A  To  A.To 

E 

Introducing  y  =  z  +  m  =  z - and  abbreviating 

A  To 

-  +  1 

k  T* 

1  k  1  1  k  1 

—  +  —  =  — ;  —  +  -  =  - we  get 

T.  A  To'  V  A  To  ( T'Y 


t  —  Of  So  —  E  f  So  —  se 

The  duration  of  the  second  period  of  movement  is  ob¬ 
tained  from  the  equation  for  z,  which  is  given  by  that 
value  of  t  for  which  2  becomes  E  once  more.  If  that 
does  not  occur  in  a  case  of  non-periodic  movement,  for 
instance,  if  the  spillway  crest  lies  below  the  level  n  —  n, 
then  the  duration  of  the  second  period  of  movement  is 
only  limited  by  a  new  occurrence  of  any  kind  of  outflow. 
Otherwise,  the  final  values  of  the  second  period  are  the 
initial  values  of  a  following  period,  which  must  be  handled 
the  same  as  the  first  case.  (Case  A.) 

The  method  of  computation  may  be  shown  best  by 
an  example.  Using  the  former  example,  we  consider  a 
spillway  of  65.7  feet  width,  the  crest  of  which  is  at  the 
static  level  n- — n.  That  is,  for  this  assumption  e 1  =  zero. 


Fig.  9. 


tPy  i  dy  y 

- 1 - +  -  =  o  85) 

It 3  T,'  dt  (T'Y 

Corresponding  to  the  investigations  regarding  the 
form  of  the  general  integral  of  this  differential  equation, 
we  must  investigate  whether  the  difference 

1  1  1 

=  - - is  positive  or  zero,  or  negative, 

(Ti'Y  (T'Y  (2  To'Y 

1  1 

which  we  obtain  by  substituting  the  values  of  —  and - 

To'  (T'Y 

1  1  k  1  k 

=  -  +  -  ( - ) 

(Ti'Y  TY  2  A  To  2  A 


The  flow  of  530  cubic  feet  per  second  is  suddenly  stopped. 
From  the  results  of  case  (A)  we  get 

zd  =  o;  td  =  106  sec.;  sd  =  +  .075  feet/sec. 

The  velocity  sd  corresponds  to  the  flow  in  the  surge  tank 

cu.ft. 

at  the  time  td  of  qd  =  .075.5380  =  404  - .  For  a 

sec. 

spillway  width  of  65.7  feet  and  for  m  =  0.6,  we  get  from 
cu.  ft. 

the  spillway  formula  q  -  =  208  h 1  dh'  and  there¬ 


fore  for 

u  .  qd  = 


sec. 
cu.  ft. 

(m  = 


282 


sec. 
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ft 


an  overfall  height  of  h'  =  1.22  feet  and  a  proportional 

sq.  ft. 

factor  k  —  3/2  .  208  .  feu*  =  345  -  and  therefore  as 

sec. 

the  distance  of  the  ideal  spillway  crest  from  the  static 
level  n  —  n  because  e1  =  zero;  E  =  .410  and  with  the 
results  of  case  a  for  z  —  E 

se  —  t-  .  073  feet  per  second 
In  order  to  determine  which  integral  formula  to  use, 
we  have 

1  1  k  1  k  1 

'TV)*  7V+  2  A  To  2  A  34.61 

1  1  1  1 

(therefore -  = - = - is  positive) 

T,2  (T1)2  (2  To1)2  (T1)2 

which  shows  that  we  should  use  the  third  form  of  the 
general  integral  which  is 


z  =  +  -354  +  i-33  e  —  l-27  e 

t  t 


165  i5-65 

s  =  —  .00805  e  +  -0810  e  which  values 

give  the  following  table : 

Seconds  t=  o  50  100  150  200 

feet  2= +.41  +1.282  +1.079  +.881  +-746 

feet 

per  sec.  5= +.073  — .00261  — .0044  — .0032  —.0024 

The  time  of  the  highest  elevation  is  determined  with 

t,  L 


*65  15.65 

s  =  o  from  equation  0  =  —  . 00805  e  +.o8ioe 

to  t,  =  40  seconds  and  2  max.  follows: 

40  40 


i65  15-65 

*  max  =  +  .354  +  1.33  e  —  1.27  e  =1.30  ft. 


t 


t 


t 


T2  T,  2  T.1 

y  =  (Ri  e  +  Ri  e  )  e 

—  E 

Further,  m  -  -  =  —  .354  feet 


A  T. 


+  1 


k  T* 


T,1  =  14.3  sec.  T1  =  51  sec.  T,  =  34.6  sec.,  so  that 

t  t  t 

+  - - - 

34.6  34.6  28.6 

2  =  +  .354  +  [Rie  +  R2e  ]e 

t  t 


This  water  level  corresponds  to  an  overflow  quantity 
of  300  cubic  feet  per  second.  For  the  determination  of 
the  ideal  overflow  height  and  the  factor  of  proportionality, 
we  used  the  maximum  overflow  quantity  of  282  cubic  feet 
per  second.  It  is  therefore  shown  that  the  assumption 
made  that  u  —  0.7  is  correct.  The  duration  of  the  over¬ 
flow  period  follows  from 

t  t 


165  15.65 

•4i  =  +  -354  +  i-33  «  —  l-27  e 

and  is  tx  =  520  seconds. 

And  therefore  the  velocity  sx  with  which  the  water  level 
reaches  the  ideal  overflow  elevation  is: 

520  520 


165  15-65 

=  +  .354  +  F,  e  +  Rz  e 


165  15.65  feet 

Si  =  —  .00805  e  +  -0810  e  =  —  .000346 - 

sec. 


t  t 


dz  R .  165  R  2  15-65 

s  =  —  = - e - e  and  for  t  =  o 

dt  165  15.65 

R,  R. 

+  .410  =  +  .354  +  f?i  +  R2;  .073  = - ; 

165  15.65 

therefore,  R,  =  1.33  feet  and  R,  =  —  1.27  feet 

t  t 


165  15.65 


With  these  initial  values  the  movement  in  the  third  period 
may  be  determined.  Because  we  assumed  that  the  spill¬ 
way  crest  is  at  the  elevation  of  the  static  level  n  —  n,  the 
overflow  period  goes  further  but  with  smaller  fluctuations 
of  water  level  and  overflow  quantities. 

If  the  spillway  is  built  not  in  the  surge  tank,  but  in 
the  main  conduit,  the  principal  equation  requires  a  sup¬ 
plement.  According  to  Fig.  10,  a  shaft  at  the  distance 
TJ  from  the  beginning  of  the  conduit  is  driven  down  to 
the  conduit  with  a  section  A1,  through  which  water  from 
the  conduit  goes  over  a  spillway  whose  ideal  overflow 
crest  is  as  before  the  distance  E  from  the  static  level 
n  —  n.  At  the  distance  L"  is  the  surge  tank  with  the 
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section  A".  In  the  hydraulic  equilibrium  with  Q,  cubic 
feet  per  second  discharge  through  the  penstocks,  the 
water  surface  in  the  overflow  shaft  will  be  below  the 
static  level  n  —  n  by  the  distance  hj  =  n'.Vi;  in  the 
surge  tank  by  hi"  ( n 1  +  n")  Vi.  Vi  is  the  velocity  in  the 
conduit  of  the  area  a  which  corresponds  to  the  discharge 
Q i.  n1  and  n"  are  the  function  coefficients  corresponding 
to  the  distances  L'  and  L" .  For  flow  fluctuations  the 
distance  Zi  and  z2  from  the  static  level  n  —  n  at  the 
velocities  v/  and  vt "  and  the  water  surfaces  in  shaft  and 
surge  tank  will  have  different  levels  and  there  will  be 
overflow  in  the  shaft  as  soon  as  the  water  surface  reaches 
the  spillway  crest.  It  is  sufficient  to  carry  through  the 
solution  of  the  principal  equation  for  the  latter  period  and 
to  handle  the  first  period  as  a  special  case  of  it. 

For  the  purpose  of  simplifying  the  work,  we  consider 
only  the  case  of  sudden  total  shut-down.  For  both  parts 
of  the  conduit,  we  get  the  following  movement  equations : 

L<  dvj 

-  +  Zi  +  n'.Vi'  =  o 

g  dt 


The  motion  equations  and  the  second  equation  for 

continuity  reach  their  values  before  as  well  as  after  the 

overflow  on  the  spillway.  The  second  equation  for  con¬ 
tinuity  is  correct  for  k  =  o,  for  periods  without  any  over- 

dzi  dz2 

flow.  If  we  consider  that  s,  =  —  and  s2  =  —  we  may, 

dt  dt 

with  the  aid  of  the  continuity  equations,  eliminate  the 
velocities  v/,  v,"  and  their  derivations.  We  get  then 
two  simultaneous  differential  equations  of  the  second 
order,  from  which  we  eliminate  again  Zi  and  its  deriva¬ 
tions,  whereas  for  the  determination  of  z2  we  get  a  linear 
differential  equation  of  the  fourth  degree  with  constant 
coefficients,  the  integration  of  which  does  not  involve  any 
great  difficulty. 

In  order  to  make  the  spillway  especially  efficient,  we 
manage  it  so  that  the  spillway  crest  lies  below  the  static 
level  n  n ,  and  this  in  a  distance  which  is  equal  to  or  a 
little  less  than  h/  =  n' .Vi.  In  this  case,  after  the  shut¬ 
down  is  finished,  a  very  rapid  overflow  will  occur  in  the 
spillway  and  will  continue,  and  finally  there  exists  a 


x 


K 


L"  .dv/ 

- h  (z2  —  Zi)  +  n"  v,”  =  0  -  (86) 

g.dt 

The  equations  which  express  continuity  follow  from 
the  consideration  that  from  the  upper  conduit  in  a  unit 
of  time  dt,  a  quantity  of  water  must  flow  into  the  shaft 
which  is  equal  to  the  algebraic  sum  of — 

i st — -the  quantity  which  flows  away  in  the  lower  con¬ 
duit,  which  is  equal  to  a.Vi" .dt. 

2nd — -the  simultaneous  filling  of  the  shaft  with  the 
quantity  A'.v/.dt. 

3rd — the  simultaneous  overflow  quantity  k  (zi — E)  dt. 
But  the  water  quantity  which  flows  through  the  lower  con¬ 
duit  in  the  time  dt  is  also  equal  to  the  simultaneous  filling 
A"  v,"  d t  in  the  surge  tank.  The  following  two  equations 
express,  therefore,  the  continuity 

a.v/  =  a.v/1  +  A’  Si  +  k  (z,  —  E ) 

and 

av/  =  AnSi  -  -  -  -  (87) 


constant  flow  over  the  spillway,  whereby  naturally  as 
much  flows  through  the  upper  conduit  as  goes  away  over 
the  spillway  and  where  the  water  level  in  the  shaft,  as 
well  as  in  the  surge  tank  lies  at  a  distance  under  n — n, 
which  corresponds  to  the  hydraulic  slope  necessary  for 
the  flow  through  the  upper  conduit. 

If  in  this  case  A'  and  the  spillway  widths  are  small 
enough  with  reference  to  the  hydraulic  slope,  to  keep  the 
preceding  fluctuations  so  near  the  static  level  of  the  water 
surface  in  the  shaft,  that  the  variation  of  the  inflow  to 
the  shaft  may  be  neglected,  the  problem  becomes  simpli¬ 
fied,  because  the  first  motion  equation  drops  out  and  in 
the  first  equation  for  continuity,  the  values  a.v/  become 
constant  and  equal  to  q. 

dzi  dz2 

Considering  that  Si  —  —  and  s2  =  — ,  we  get  from 
dt  dt 

the  first  equation  of  continuity  (87) 


INFLUENCE 


OF  h  =  u.  v 2. 
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dzz  dsi 

A"  . -  -r  A> -  +  k  (zi  +  E)  —  q.  (88) 

dt  dt 

and  the  second  motion  equation  becomes  zero. 

L"  d2z2  dz2 

—  A" - r 2  A"  -  +  a  (z2  —  Zi)  =  o  (89) 

g  dt 2  dt 

From  this  last  equation,  Zi  and  its  derivative  may  be 
determined  and  substituted  in  equation  30c.  For  the 
solution  of  z2  we  get  a  linear  differential  equation  of  the 
third  degree  which  may  readily  be  integrated.  The  re¬ 
sults  of  such  an  investigation  are  apparent  in  Fig.  1 1  ;  we 
see  that  at  the  time  of  the  highest  water  elevation  in  the 
surge  tank  the  height  of  water  flowing  over  the  spillway 
is  directly  proportional  to  the  height  of  the  water  in  the 
surge  tank  above  the  static  level,  but  the  quantity  of 
water  is  greater  during  the  time  the  level  in  the  surge 
tank  is  falling  and  continues  until  the  level  in  the  surge 
tank  has  reached  its  lowest  elevation.  The  extent  of  the 
fluctuations  is  not  only  dependent  upon  the  conduit,  surge 
tank  and  spillway  dimensions,  but,  naturally,  also  de¬ 
pends  upon  the  area  A'.  The  complete  computation  of 
the  case  would  exceed  the  space  available  here. 

Investigation  of  the  Influence  of  h  =  u.v2. 

Finally,  a  study  may  be  added  which  permits  a  com¬ 
parison  between  the  results  of  the  assumption  already 
made  (that  h  —  n.v),  and  the  results  obtained  if  we  use 

Lv2 

the  usual  formula  h  =  u.v2  =  u1 . - where  the  friction 

P.2g 

coefficient  u1  is  again  assumed  as  constant  (P  =  hydraulic 
radius  of  the  conduit).  With  this  condition  it  is  neces¬ 
sary  to  use  the  first  original  equation  (15)  in  which  we 
must,  according  to  the  direction  of  the  flow  in  the  con¬ 
duit,  be  careful  to  use  the  right  sign  for  h.  If  this  makes 
the  investigation  somewhat  complicated,  it  may  be  said 
that  the  differential  equation,  which  we  obtain  introducing 
the  value  of  v  in  the  continuity  equation  (17)  is  not  only 
of  the  second  (i.e.,  higher)  order,  but  also  of  the  second 
degree,  so  that  difficulties  occur  in  the  determination  of 
the  general  integral.  This  is  especially  the  case  if  in  the 
equation  mentioned  q  >  o  and  variable. 

Therefore,  in  the  following  we  consider  only  for  com¬ 
parative  purposes  the  case  of  a  sudden  shut-down  of  Q 
and  for  the  phase  of  the  first  surge.  In  order  to  get 
more  generality,  we  first  consider  the  surge  tank  section 
as  variable  with  the  height,  which  means  that  it  is  as¬ 
sumed  as  a  given  function  of  z.  The  motion  equation  is 
then  : 

L  dv 

- h  z  +  h  =  o.  -  -  (15  A). 

g  dt 

The  equation  of  continuity  is  : 

a.v  =  A  .s  -  -  -  (17  B) 


function  of  a  variable  x  where  z  =  x  —  hi  and  where  V  is 
measured  positive  up  from  that  level  which  lies  below 
n  —  n.  That  is,  for  z  =  —  hi  or  x  =  o ;  V  —  o. 
u 1  Lv2  u1  L  v2 

With  h  =  -  and  hi  =  -  and  multiplying  (90) 

P.2g  P . 2g 


2g 

with - we  get 

L  a 

2  g  ul 

d(vy  +  [—  .  x  +  —  (v2  —  V22)}  dV  =  o 
L.a  P.a 

2  g  u‘  1 

If  we  say  v2  =  y  v2*  =  y2  -  =  a,  —  =  —  the 

L.a  P.a  a2 

total  differential  equation  follows : 

1 

dy  +  [a,  x  +  —  (y  —  y2)j  dV  =  o  which  integration  can 
&2 


V 


be  done  introducing  the  integrating 
V  V 


the  values  for  y,  a,,  a2. 


&2 

factor  m  =  e 


and  introducing 


ulV 


2  g  P.a 

v2  =  V21 - e 

L.a 


u1  V 
V  +  - 

P.a 

x  e  .  dV  (91) 


With  help  of  this  equation  v  may  be  found  as  a  function 
of  V  and  x,  because  by  given  relation  of  the  volume 
V  to  x  and  vice  versa,  in  formula  (91) 


ul  V 


P.a 

.dV  may  be  determined  by  means  of 


quadrature  even  if  V  =  F  (x)  is  given  only  by  a  curve. 
The  method  is  illustrated  by  the  diagram  shown  in  Fig.  12. 
The  relation  of  the  volume  V  to  x  is  shown  by  a  curve, 
the  abscissae  of  which  represent  the  values  of  V  and  the 
ordinates  the  values  of  x.  The  construction  of  the  X  curve 

I  —  I  may  be  accomplished,  based  on  the  given  dimensions 
of  the  surge  tank.  V  =  o  corresponds,  as  said  before, 
to  the  operating  level  for  a  flow  of  Q 1  cubic  feet  per 
second.  With  the  same  abscissae,  we  construct  a  curve 

II  —  II.  the  ordinates  of  which  are  given  by  the  values 

ul  V 

+  - 

P.a 

e  ■  and  another  curve  III  —  III  with  the  ordinates 

ulV 


If  we  multiply  equation  15  A  with  a.v.dt,  it  follows,  con¬ 
sidering  equation  (17  B) 

L.a 

- v.dv  +  z  A  s  dt  +  h.A.s.dt  =  0  (90) 

g 

A  .s  dt  =  A  .dz  =  dV  represents  the  change  of  the  water 
volume  V  in  the  time  dt.  It  also  represents,  with  given 
surge  tank  dimensions  as  function  of  z  or  with  reference 
to  the  graphical  demonstration  appearing  further  on,  the 


P.a 

x.e  and  then  the  integral  curve  IV  —  IV.  If  we  con¬ 
struct  a  further  curve  V  —  V,  the  ordinates  of  which  we 
get  if  we  multiply  the  ordinate  values  of  the  integral  curve 
u1  V 


2  g  P.a 

by  - e  and  if  we  construct  at  the  distance  (vt*)  a 

L.a 
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line  parallel  to  the  axis  of  abscissae  we  then  see  that 
the  difference  of  the  ordinates  between  the  (v22)  line  and 
the  curve  V  —  V  gives  the  values  of  v 2.  Naturally,  we 
must  select  suitable  scales  for  all  curves,  that  is  to  say, 
for  their  ordinates.  The  scale  for  the  last  curve  is  natur¬ 
ally  the  same  as  for  (7;22).  The  intersection  of  the  last 
curve  with  the  parallel  (i/22)  determines  the  ordinate  which 
gives  the  value  of  V  on  the  axis  of  abscissae,  for  which 
v  =  o,  that  is,  when  the  highest  water  elevation  occurs. 
T  he  corresponding  ordinate  of  the  x:  curve  gives  the 
height  of  the  highest  elevation  above  the  operating  level. 

If  A  is  constant:  V  =  A.  x  and  equation  91  becomes 
A.x  Ax 


—  m1 -  x  u' 


2  g  A  1  A  1 

If  we  introduce  -  .  —  =  — ;  u1  -  =  —  we  get 

L  a  B2  P.a  D 


1  2  g  A  1 


— 

-  — 

B  = 

B2 

L.a 

2.094 

v22  B2 

6.632  x  1.447s 

+  1 

1  ( - 

-)  + 

D2  9.6s 

Therefore 

X 

X  D 

—  +  e  -  -  2.00 

D 

X 

—  =  1.86  X  —  17.81  feet 

D 

In  case  (A),  where  we  assume  the  same  conditions 
with  hi  =  n.v  in  which  n  =  1.445  seconds,  we  got 
z  max  =  6.55'  and  therefore  the  highest  elevation  above 
the  initial  level  z  max  +  hi  =  6.55  +  9.6  =  16. 15'.  This 


value  is  smaller  than  the  value  which  we  computed  for 
X  by  1.66  feet,  and  it  is  necessary  for  us  to  correct  the 
formulae  obtained  under  the  assumption  that  h  =  n.v,  in 
order  to  get  as  near  as  possible  to  the  new  results.  For 
this  purpose,  we  assume  that  for  the  expression  of  the 
elevation  x  above  the  initial  level  the  formula 


D2  x  D 

v2  =  v  2  +  -  (1 - e  )  (92) 

B2  D_ 

And  for  the  computation  of  the  maximum  value  X 

of  x  (v  =  o) 

X 

vA.B2  X  D 

-  +  1  =  —  +  e  (93) 

D2  D 

The  results  here  obtained  may  now  be  used  for  com¬ 
parative  purposes  with  the  results  of  the  example  of 

case  (A). 

v2 

v 2  =  6.63  feet  per  sec.  -  =  .68  feet 

2 

w1 

hi  —  9.6  feet  =  —  .  9050.0.68 

P 


t 


2  To  t  ' 

x  =  h,  +  z  =  hi  +  R.e  sin  f/3  4 - ) 

Ti 


t 


(94) 


dx  R 

dt  T 


2  T„  t 

sin  (27 — 0 - ) 

T, 


and  therefore 

t 


d2x  R 

dt2  T2 


2  To  t 

sin  (27  —  0 - ) 

T, 


2  To 

with  tg  7  = - 

T, 

and 

11  1 

T2  T2+  (2  To)2 

(95) 


should  be  used,  where  the  values  R  and  0  are  again  to 
be  determined  from  the  initial  conditions,  that  is  to  say, 


u* 


1 


r  u'  A 


P  641 


D  Pa  9.6 


dx 

where  t  =  o,  x  =  o  and  z  =  ht,  —  =  c,,  but  the  values 


INFLUENCE  OF  h  =  u.v\ 


2  3 


♦ 


T„  and  T  and  T,  must  be  determined,  not  as  before  from 
the  given  dimension,  but  from  the  condition  that  X  be¬ 
comes  the  maximum  value  of  the  elevation  above  the 

dx 

initial  elevation,  that  is,  for  x  =  X  the  derivative  —  must 

dt 

become  zero,  and  as  before  for  x  =  o,  the  second  de- 
d2x 

rivative  -  becomes  zero.  From  the  latter  condition 

dt 2 

follows  : 

27  —  &  =  i8o°  =  Jr  therefore  7  —  B  =  180  — -7 
B  =  27  —  i8o°  and  for  t  =  o  we  get 

o  =  hx  —  R  .  sin  27 


because  —  = -  a!  ( - )2  +  1.  See  equation  (38). 

T  2  To  V  Tx 

dx 

The  time  f*  when  x  =  xmu  —  X  ( —  1  o)  is  computed  from 

dt 

tx 

—  =  7  —  &  =  ir  —  7.  Therefore 
Tx 

T 1  tx 


2  To 

X  =  hx  +  Re 


Tx 

and,  as -  =  ctg  7 

2  T# 


X  —  hx  +  R  e 


Tx 

sin  (B  +  7  —  0) 


—  —  7)  ctg  7 

.  sin  7 


(96) 


cases,  must  have  the  same  values.  If  these  values  are 
not  exactly  equal,  but  if  we  get  approximately  equal 
values  it  depends  on  the  degree  of  approximation  desired 
whether  we  use  the  supplementary  formula  (92)  or  not. 

In  order  to  investigate  this  case,  we  proceed  to  com¬ 
pute  the  example. 

With  X  =  17.81  feet  hx  —  9.6  feet  we  get 

+  («■  —  7)  cotg  7  1 

cos  7  e  -  -  =  .581  7  =  710  50' 

X 

2( - 1) 

hx 

2  T. 

0  =  27  —  1800  -  —  36°  20' ;  tgy  —  -  —  3.047  ; 

Tx 

hi  feet 

R  =  -  =  16.2  feet  with  c,  =  0.098  - 

sin  2  7  sec. 


2  To 


R 

—  sin  7 

Ci 


R 

tg2  7  +  1  =  —  tgy 

Cx 


300  sec. 


T,  =  164.3  sec-  T  -  155.9  sec.  tz  (x  =  X)  =  310  sec. 

t 


500  t 

and  therefore  x  =  9.6  +  16.2  e  sin  ( - 36°) 

164.3 

t 


16.2  500  t  s.i4  5380 

5  = - e  sin  (71 0  50'  + - );v  = - =  sx - 

r55-9  164.3  a  80 


and  we  can  first  determine  7  with 

hi  +  («■  —  7)  cotg  7 


With  these  formulae,  we  may  compute  the  relative 
values  of  x,  s  and  v  for  the  different  values  of  t.  From 
the  supplementary  formula  (92) 


R  =  - ;  cos  7  .  e  =  - 

sin  27  X 

2  ( - 1) 

hi 


Is  7  with  this  equation  computed,  we  have  for  the 
other  values 


These  values  may  then  be  introduced  in  the  formulas  for 
dx 

x  and  —  =  s  -  -  -  (94) 

dt 

If  the  movement  so  determined  must  be  similar  to 

L  v 2 

that  which  occurs  if  h  =  u'  —  .  — ,  then  the  velocities 

P  2  g 
A 

v  =  s  — ,  which  we  get  for  the  same  elevation  x  in  both 

a 


x 

D*  x  D 

v2  =  vS  +  —  (1 - e  )  .  (o2) 

B2  D 

we  may  compute  the  values  of  v  for  the  same  values  of  x 
or  they  may  be  determined  graphically  from  Fig.  12. 

L  v2 

They  correspond  to  the  assumption  that  h  =  u'  _ .  — ■ 

P  2  g 

and  we  get  the  following  table  of  values : 


t  Sec. 

0 

50 

100 

ISO 

200 

250 

300 

310 

x  feet 

0 

4.06 

9.21 

12,9 

15.7 

17.4 

17.8 

17.81 

s  feet  , 

;  sec. 

098 

.092 

.082 

.0*23 

.0425 

.023 

.00328 

.00 

and  from 

that 

v  feet  I 

'  sec. 

6.63 

6.25 

5.49 

4.23 

2.92 

1.51 

.20 

.00 

But  from  Formula  (92) 

v  feet  , 

'  sec. 

6.63* 

6.3 

5.38 

4.05 

2.83 

1.35 

.066 

.00 

If  we  plot  in  the  rectangular  co-ordinate  system  (Fig. 
12)  the  two  values  of  v  as  ordinates  with  the  values  of  x 
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as  corresponding  abscissae,  then  we  see  from  the  curves 
thus  obtained  that  we  may  use  the  supplementary  formula 
with  sufficient  accuracy  for  the  determination  of  the  move¬ 
ment.  If  we  compare  the  values  of  To  and  T  as  they  have 
been  found  for  both  assumptions  of  h,  that  is  to  say,  for 
To  the  values  194.5  and  250.0  and  for  T  the  values  137.5 
and  155.9,  then  we  see  that  the  second  value  of  T0  is  1.28 
times  as  large  as  the  first,  the  second  value  of  T  is 
1. 13  =  vG.28  times  as  large  as  its  first  value.  The  first 

L 

values  have  been,  found  from  the  formulae  To  =  - and 

_  n.g 

IT  A 

T=  4  /  —  .  —  and  from  this  we  see  the  importance  of  the 

\  g  a 

supplementary  formula  (92).  It  is  evolved  from  the  principal 
equation  15  if  we  introduce  instead  of  L,  a  distance  L* 
which  is  28%  larger  than  L  (provided  that  the  value  n  is 
the  same  according  to  the  relation  lu  =  n.v,)  then  with 
the  distance  L* ,  we  compute  the  values  T  and  T0.  It  is 
evident  that  with  the  introduction  of  the  formula  of  re¬ 
sistance  h  1  =  n.v,  the  influence  of  the  friction  has  been 
introduced  in  the  equation  as  too  large.  We  have,  there¬ 
fore,  assumed  in  the  supplementary  formula  the  working 
capacity  and,  therefore,  also  the  mass  of  the  conduit 
volume,  larger  than  they  really  are,  but  obtained  the  right 
value  of  the  elevation  through  the  supplementary  formula. 
The  friction  height  which  corresponds  to  the  elevation 
X  =  17.81  feet  can  be  obtained  by  computing  the  work 
balance,  according  to  the  former  scheme,  if  we  introduce 
the  actual  length  of  the  conduit. 

G  becomes  5,380  x  17.81  =  2,994  tons; 

yi  =  2  17.81  =  8.9  feet. 

USEFUL  WORK  IN  FEET  TONS  USED  WORK 

15,800  KINETIC  ENERGY  OF  CONDUIT  VOLUME 

28,740  POTENTIAL  EN ERGY  OF  G  =  2994t  FOR  h  =  9-6 

LIFTING  WORK  OF  2990t  FOR  8  9'  26,650  t  Total 

Total  44,540  FRICTION  WORK  17,890/  44,510 

The  total  friction  work  of  17,900  foot-tons  corresponds, 
therefore,  to  an  average  friction  height  with  respect  to  the 
weight  of  the  water  of 

17,900 

H  average  =  -  =  .  60  h, 

2,994 

For  case  (A)  we  found  h  average  =  .755  h,  as  the 
average  value  which  corresponds  to  an  elevation  of  only 
16.75  feet  actually.  The  maximum  elevation  above  the 
initial  level  will  be  larger  than  16.75  ^eet  and  smaller  than 
17.81  feet,  therefore,  0.60  h,  <  h  average  <  .755  h,.  The 
true  value  will  be  apparent  if  observation  results  in  con¬ 
structions  of  larger  dimensions  than  can  be  obtained. 
We  may  remark  also  that  the  above  determined  values 
of  T  and  T0  in  the  supplementary  formula  are  only  correct 
for  the  first  surge  from  x  =  o  to  x  —  X.  We  may  carry 
the  mathematical  investigation  further  and  find  the  values 
for  the  drop  from  the  first  maximum  elevation  to  the  first 
minimum  level.  Other  values  for  T  and  T»  must  be  intro¬ 
duced,  which  correspond  to  smaller  values  of  L*.  This 
repeats  itself  for  the  further  phases  of  increase  and  de¬ 
crease  until  finally  L*  becomes  nearly  equal  to  L.  This 
corresponds  to  the  decrease  of  the  velocity  fluctuations 
during  the  movement  and  also  with  the  decrease  of  the 
influence  of  friction. 

Conclusions. 

We  may  summarize  our  results  as  follows : 

1.  If  the  conduit  ends  in  a  pond,  from  which  the  pen¬ 
stock  goes  directly  to  the  turbine,  and  if  the  area  of  the 


pond  surface  A  >  150.771. a  where  a  —  the  conduit  area  in 
square  feet  and  m  the  number  of  miles  conduit  length  : 
we  need  not  expect  any  periodical  fluctuations  of  the  out¬ 
flow,  even  if  the  outflow  varies  during  the  observed  time. 
If  A  <  150  ma  as  is  the  case  with  artificial  surge  tanks, 
such  periodical  fluctuations  may  occur. 

2.  The  maximum  elevation  above  the  initial  level 
(operation  level  for  full  outflow)  has  approximately  the 
same  ratio  for  a  short  time  for  a  sudden  as  well  as  for  a 
gradual  close  of  the  whole  outflow. 

3.  The  value  of  the  maximum  elevation  above  the 
initial  level  may  best  be  computed  by  means  of  the  work 
balance,  where  for  the  average  friction  height  0.7  of  the 
value  of  that  height  h ,  must  be  introduced,  which  corre¬ 
sponds  to  the  drop  of  the  water  surface  in  the  surge  tank 
below  the  level  of  the  water  surface  before  the  conduit  for 
full  outflow. 

4.  The  extent  of  drop  below  the  static  level  for  full 
opening  is  approximately  equal  to  the  rise  above  computed. 

5.  For  the  determination  of  the  oscillations  during 
and  after  completed  shut-down,  we  may  use  for  a  fair 
degree  of  approximation  the  formulae,  computations  and 
demonstration  methods  which  correspond  to  the  analysis 
above  developed  on  the  theory  of  damped  oscillations. 
The  values  of  T  and  T0  and  therefore  also  T,  may  be  de¬ 
termined  if  investigating  the  maximum  elevations,  as 
shown  in  our  last  study. 

6.  For  operation,  which  results  in  a  periodical  varia¬ 
tion  of  the  outflow,  resonance  features  may  occur. 

7.  The  maximum  elevations  may.  be  decreased  by 
construction  of  spillways  in  the  surge  tank  or  in  the 
conduit. 

The  theory  furnishes  a  good  explanation  concerning 
the  manner  and  the  relative  value  of  the  operations  close 
enough  for  most  practical  purposes.  But  a  high  degree 
of  accuracy  will  be  accomplished  only  where  experiments 
furnish  much  more  exact  values  for  friction  in  conduits 
under  different  conditions. 

Appendix  A — Surge  Tank  With  Variable  Cross=section. 

Assuming  a  quadratical  relation  between  friction- 
head  and  conduit  velocity  the  following  relation  was  ob¬ 
tained.  (See  Formula  91.) 


u1  V  u'V 


The  maximum  surge  is  obtained  by  making  v2  =  o 
and  solving  the  latter  equation  for  X.  The  previously 
described  graphical  method  may  be  applied  without 
further  comment  to  surge  tanks  of  variable  cross-section 
with  the  exception  that  the  X  curve  is  no  longer  a  linear 
function  of  V.  For  this  reason  and  because  of  complex 
integration  due  to  the  insertion  of  a  complicated  function 
the  graphical  method  is  more  appropriate  than  the 
analytical. 

In  the  following  example  this  method  is  applied  to  a 
surge  tank,  circular  in  section,  which  has  a  sectional  area 
of  5,380  sq.  ft.  at  the  operating  level,  as  was  assumed 
before.  However,  in  this  case  the  sectional  area  in¬ 
creases  in  direct  proportion  to  the  height  in  such  a  ratio 
that  at  a  height  of  19.7  ft.  the  sectional  area  is  6,450  sq. 
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ft.  (See  Fig.  14.)  All  other  data  for  the  problem  is  the 
same  as  assumed  in  the  main  article. 

L  =  9.°5°  ft- 
a  =  80  sq.  ft. 
p  -  32.8  ft. 

Q  =  ^30  cu.  ft.  per  sec. 
v2  =  6.63  ft.  per  sec. 

P  =  2.44 
hi  =  9.6  ft. 
u'  —  .0038 


P.a 


5I>3°° 


CL2 


64.4 


51 >3°° 


V 

V  + - 

5G300 


x  e 


dV 


6.63*  =  44.0 

9050.80  J  0 

With  this  graphical  method  (Fig.  13)  a  maximum 
surge  of  only  17.0  ft.  is  obtained,  whereas  in  the  former 
example,  with  a  surge  tank  of  a  constant  sectional  area 


Fig.  14. 


Fig.  15. 


Fig.  Ifi. 
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of  5,380  sq.  ft.  a  surge  of  17.8  ft.  resulted.  '  But  the 
volume  of  water  lifted  is,  in  this  case,  somewhat  greater, 
that  is,  98,800  cu.  ft.,  as  compared  to  but  95,600  cu.  ft. 
in  the  former  case. 

In  the  same  easy  way  the  graphical  method  may  be 
applied  to  a  surge  tank  in  which  the  cross-section  varies 
irregularly  with  the  height.  For  instance,  as  the  Albula 
plant  of  the  city  of  Zurich.  (Fig.  15.)  In  this  a  break 
occurs  at  the  points  on  the  %  curve  where  the  ratio  be¬ 
tween  sectional  area  and  height  changes.  (See  Fig.  16.) 

Appendix  B — Simplified  and  Approximate  Formula  for  the 
Determination  of  Surge  Tanks  With  Constant 
Cross=Section. 

By  A.  Strickler. 

The  designing  engineer  must  often  consider  several 
solutions  for  one  problem  and  from  these  select  the  best 
for  the  actual  work.  For  these  preliminary  investigations 
he  is  in  need  of  a  simple  slide-rule  formula.  In  the  follow¬ 
ing  such  a  formula  will  be  derived  for  the  required  area 
of  a  surge  tank  with  given  maximum  surge.. 

As  initial  formula,  equation  51  may  be  used  for  a  sud- 


M  .  w,1 

+  G  .  h,  =  G  .  y,  +  A1  -  (51) 

2 

The  notation  being  the  same  as  before,  the  average 
friction  head  in  the  surge  tank  for  the  first  surge  may  be 
called  ha,.  Equation  51  can  now  be  written: 

-CL.  L  y  max 

•  v'2  +  G  ( - h,)  +  G  .  ha  =  o 

2  g  2 


is  a  regular  cosine  curve.  This  would  be  true  if  the  water 
level  fluctuations  were  an  undamped  harmonic  of  the  form 
t 

z  —  C  .  sin  ( — )  (see  Fig.  17). 

T 

Making  the  above  assumption 


5 .  A  A  dz  A  C 

t 

v  =  -  =  —  .  —  =  —  .  — 

.  cos  ( — ) 

a  a  dt  0  T 

T 

AC  0  A 

C 

and  also  Vi  =  —  .  —  .  cos  —  =  — 

.  - 

a  t  t  a 

T 

t 

t 

v  —  Vi  cos  — ;  h  =  hi  .  cosa  ( — ) 

T 

T 

Substituting  these  values  in  equation  (51b)  there  results: 

«■ 

—  T 

r  1  u  t  hi 

C 2  t 

/  ha=  —  K  .  cos1  ( — )  dt  — -  J 

cos1  ( — )  dt  ( 51c) 

J  ti  0  T  n  J 

- T 

2 

t  T 

hi  n  .  t  hi 

hi 

=  -  x  -  =  — 

ha=  — 

JT  4  2 

- .T 

2 

The  average  friction  head  in  the  tank  for  the  time  up 
to  the  maximum  surge  is  therefore  equal,  under  the  above 
mentioned  assumption,  to  half  the  friction  head  which 
exists  in  the  conduit  during  operation.  This  value  of  ha 
introduced  into  equation  51a  (work  balance)  gives: 
a.L  A  A.  ymax. hi 

- .  Vi*  +  -  .  yJmax  -  A  .  y max  .  hi  +  - =  0 

2  g  2  2 

L  L  .  V* 

—  .  v*  a - 

g  S 


or  A  =  a  . 


ymax  (ymax  -  hi) 

a.v *  2  L 


(z max  -+■  hi)  Zmax 


.-.  A  = 


2  g  (z max  +  hi)  Zm ax. 


This  formula  gives  results  somewhat  too  large  (i% 
to  3%)  but  is  exact  enough  for  preliminary  computations. 
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t,  =  time  from  the  beginning  of  the  surge  to  the  moment 
of  maximum  elevation.  In  order  to  integrate  equation 
51b,  v  must  be  known  as  function  of  t.  The  simplifying 
assumption  is  now  made  that  the  first  part  of  the  v  curve 
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